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Abstract Observations of thermal radiation from neutron stars can potentially pro¬ 
vide information about the states of supranuclear matter in the interiors of these stars 
with the aid of the theory of neutron-star thermal evolution. We review the basics of 
this theory for isolated neutron stars with strong magnetic fields, including most rel¬ 
evant thermodynamic and kinetic properties in the stellar core, crust, and blanketing 
envelopes. 
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1 Introduction 


The first works on neutron star cooling and thermal emission (Stabler 

196(4 Tsurutal 

1964: Chiu and SalDeteJl964: Mortonlll964: 

Bahcall and Wolflll965al 

g) appeared at 


the epoch of the discoveri es of X-ray sources outside the Soly Syste m in the rocket 
and balloon experiments ( Giacconi et al. 1962 : Bowver et al. 1964allb ). The authors 


estimated cooling rates and surface temperatures in order to answer the question, 
whether a neutron star can be detected in this kind of experiments. However, the 
first attempts failed to prove t he relation between ne utron stars and newly discovered 


X-ray sources. In particular, iBowver et al.l (Il964bl) measured the size of the X-ray 


source in the Crab Nebula from observations during a lunar occultation on July 7, 
1964. Their result, ^ 10*^ km, indicated that the source was much larger than a neu¬ 
tron star should be. Ironically, there was a neutron star th ere, the famous Cr ab pulsar, 
but it was hid den within a compact plerion pulsar nebula. Kardashev ( 1964ll and later 


Pacinil (1 19671) conjectured that the Crab Nebula could be powered by the neutron-star 


rotational energy, which was transferred to the nebula via the magnetic field, but this 
model remained a hypothesis . Curiously, the Crab puls ar was observed as a scintillat¬ 
ing radio source s ince 1962 (Hewish and Okovjl 196.5 ), but the nature of this source 


remained unclear. Sandage et al. ( 19661) identified Sco X-1, the first detec ted and the 


bright est cosmic X-ray source, as an optical object of 13th magnitude. IShklovskv 


( 19671) analyzed these observations and concluded that the X-ray radiation of Sco 


X-1 originated from the accretion of matter o nto a neutron star from its com panion. 


Later this conjecture was proved to be true (Ide Freitas Pacheco et al.l 1 19771) . but at 
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the time it was refuted ( Camerorilll967 ). Because of these early confusions, the first 
generall y accepted evidence of neutron stars was provided only by the discovery of 


pulsars (IHewish et al .11 19681) after a successf ul competiti on of the theoretical inter¬ 


pretation of pulsars as rotating neu tron stars ( Gold ll968h with numerous alternative 


hypotheses (see, e.g., the review by Ginzbur j 


19711). 


The foundation of the ri gorous cooling theory was laid by Tsurutal ( 19641) and 
Tsuruta and CameronI (119661) . who formulated the main elements of the theory; the 
relation between the internal and surface temperatures of a neutron star, the neutrino 
and photon cooling stages, etc. After the discovery of neutron stars, a search for their 
soft X-ray thermal emission has become a topical challenge, which stimulated the 
de velopment of th e co oling theory. The first decade o f this development was reviewed 
by (122^ and Nomoto and Tsurutal ( 1981 ah . 


Thornel ( 1977 ) presented the complete set of equations describing the mechani¬ 


cal and thermal structure and evolution of a spherically symmetric star at hydrostatic 
equilibrium in the framework of General Relativity (GR). The GR effects on the 
the rmal evolution of neutron s t ars were first included into th e cooling calculations 
by Glen and Sutherland! ( 198^ : Nomoto and Tsurut3 ( 1981b ): Van Riper and Lam3 
(Il98lh . A generally relativist ic cooling code for a spherically symmetri c non- 
barot ropic star was written by Richarson et al. ( 1982h . Nomoto and Tsuruta ( 1986 . 
1987h studied neutron star cooling using an updated physics input and discussed the 
role o f different phys ical factors for thermal relaxation of different models of neutron 


stars. ITsuruta ( 19861) provided a comprehensive review of the neutron star cooling 
theory with a comparison of the results of different research groups obtained by mid- 
1980s. 

The early studies of the neutron-star cooling were mostly focused on the stan¬ 
dard scenario where the neutrino emission from the stellar core was produced mainly 
by the modified Urea (Murca) processes, which compete with neutrino emission via 
plasmon decay, nucleon bremsstrahlung, etc. The enhanced (accelerated) cooling due 
to the direct Urea (Durca) processes was belie ved possible only if the core con¬ 
tains a pion condensate or a quark plasma (e.g., Tsurutal 19791 : Glen and Sutherland 
1980l: Van Riper and Lamblll981 ; Richarson et al.lll982l). Bv the end of 1980s a new 


coolin g agent, kaon condensate, was introduced ( Brown et aDll988 ; Page and BaronI 

1990l). The studie s of th e enhanced cooling were intensified after the discovery 


by iLattimer et al.l (119911) that the Durca process is allowed in the neutron star 


core with the standard nuclear composition for some realistic equations of state 
(EoS) without “exotic” models. The standard and Durca-enhanced neutron star cool- 
ing scenarios were compa red in a number of numerical simulations starting from 
Page and Applegate ( 19921) . who also noticed that nucleon superfluidity becomes the 
strongest cooling regulator in the Durca-allowed stellar kernels. This result triggered 
a flow of papers on the cooling of superfluid neutron stars. 

The progress in the theoretical studies of the neutron-star thermal evolution was 
influenced in the 1980s and 1990s by the spectacular progress of the X-ray astronomy, 
notably due to the space observatories Einstein (1978-1981), EXOSAT (1983-1986), 
and ROSAT (1990-1998). ROSAT was the first to reliably detect X-ray thermal ra¬ 
diation fro m isolated neutr o n stars. This theoretical and observationa l progr ess was 
reviewed by Tsuruta ( 1998 ): Yakovlev and Pethick (2004); Page et al. (12004 ). 
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In the 21st century, the data collected by X-ray observatories Chandra and XMM- 
Newton give a new impetus to the development of the cooling theory. Some new 
theoretical results on the coo ling of neutron stars a n d rela ti on of t he theory to ob - 
servations were reviewed by Yakovlev et al.l ( 20081) : Page (2009); Tsuruta ( 2009 ). 
Recently, 2D simula tions of the fully coupled thermal an d magnetic field evolution 
have been possible ( Pons et ^ 2009; Vigano et al.l 2013 ). mostly motivated by the 
increasing number of observed magnetars and high magnetic field pulsars. 

The theory of thermal evolution of neutron stars has different aspects associ¬ 
ated with rotation, accretion, etc. In this review, we will mostly focus on the physics 
that determines thermal structure and evolution of slowly rotating non-accreting neu¬ 
tron stars, whose thermal emission can be substantially affected by strong magnetic 
fields. The processes of formation of thermal spectra in the outermost layers of such 
stars are explicitly excluded from this p aper but considered in the companion re¬ 
view (IPotekhin, De Luca, and Ponsll2QL5l hereafter Paper I). We will pay a special 
attention to the effects of strong magnetic fields on the thermal structure and heat 
conduction in the crust and heat-blanketing envelopes of neutron stars. 


2 The essential physics of neutron star cooling 

In this section we briefly present the essential physical ingredients needed to build 
a model of a cooling neutron star regardless its magnetic field. The effects of strong 
magnetic fields will be discussed in subsequent sections, starting from Sect. |4] 


2.1 Structure and composition of a neutron star 

A neutron star is born hot (Ri lO" K) and lepton-rich, but only a few days after its 
birth, its temperature drops to a few x 10^ K. Thus, the Fermi energy gp of all particles 
is much higher than the kinetic thermal energy in most of the star volume, except in 
the thin outermost layers (a few meters thick), which does not affect the mechanical 
and thermal structure of the rest of the star. Therefore, a good approximation is to 
describe the state of matter as cold nuclear matter in beta equilibrium, resulting in an 
effectively barotropic EoS. The mechanical structure of the star is decoupled from its 
thermal structure and can be calculated only once and kept fixed during the thermal 
evolution simulations. 

To a very good approximation, the mechanical structure can be assumed to be 
spherical. Appreciable deviations from the spherical symmetry can be caused by 
ultra-strong magnetic fields (Z? > 10^^ G) or by rotation with ultra-short periods (less 
than a few milliseconds), but we will not consider suc h extreme cases. T hen the 
space-time is described by the Schwarzschild metric (e.g. Misner et al.lll973l) 


d.v^ = —e^‘*’(''^c^df^-|-e^^^'')d4-(-r^(d0^-|-sin^ 


( 1 ) 


where (r, 0,^) are the standard spherical coordinates, = 1 — 2GMr/c^r, and 

0(r) is determined by equation 


d0(r)/dP(r) = - [P{r) + p{ry] 


,21 -1 


( 2 ) 
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with the boundary condition = 1 — rg/R at the stellar radius R. Here, Kg = 

2GM/c^ = 2.95 (M/Mq) km is the Schwarzschild radius, M = Mr is the stellar mass, 
Mr = 4 k Jq p{r)r^dr is the mass inside a sphere of radius r, G is the gravitational 
constant, c is the speed of light, P is the pressure, and p is the gravitational mass 
density. 

The mechanical structure of a spherically symmetric star is described by the 
Tolman-Oppenheimer-Volkoff equation 


dr r^ \ Pc^/ V MrC^ J \ rc^ J 


(3) 


where r is the radial coordinate measured from the stellar center. In order to determine 
the stellar mechanical structure, Eq. Q should be supplemented by an EoS, which 
depends on a microscopic physical model (Sect. Several qualitatively different 
regions can be distinguished in a neutron star, from the center to the surface; the inner 
and outer core, the mantl e, the inner and outer crust, the ocean, and the atmosphere 
(e.g., Haensel et alJl2007l) . 

The outer core of a neutron star has mass density 0.5po <P< 2po, where po = 
2.8 X 10^^ g cm^^ is the nuclear saturation density (the typical density of a heavy 
atomic nucleus). It is usually several kilometers thick and contains most of the stellar 
mass. The outer core is mostly composed of neutrons with an admixture of the protons 
and leptons - electrons and muons (npep matter). 

The inner core, which can exist in rather massive neutron stars, M > 1.5 Mq, occu¬ 
pies the central part with p > 2po. It is defined as the region where the composition is 
uncertain, but probably more rich than simply neutrons and protons. Its composition 
and properties are not well known because the results of their calculation strongly de¬ 
pend on details on the theoretical model of collective fundamental interactions. Some 
of the proposed models envision the following hypothetical options; 


1. hyperonization of matter - the appearance of various hyperons (first of all. A- and 

-hypemns-npepAE matter); 

2. pion condensation - formation of a Bose condensate of collective interactions 
with the properties of ;r-mesons; 

3. kaon condensation - formation of a similar condensate of /T-mesons; 

4. deconfinement - phase transition to quark matter. 

The last three options are often called exotic 
paper we will not consider the exotic matter in any detail. 

In the stellar crust and ocean the matter is less extraordinary; it contains elec¬ 
trons, nuclei composed of protons and neutrons, and, in the inner crust, quasi-free 
neutrons. Nevertheless, this region is also under extreme conditions (density, temper¬ 
ature, magnetic field) that cannot be reproduced in the laboratory. In the crust, which 
is normally ^ 1 km thick, the nuclei are arranged into a crystalline lattice, and in the 
ocean with a typical depth from a few to ^ 100 meters (depending on temperature) 
they form a liquid (see Sect. 12.4.1b . 

With increasing density, nuclei become progressively neutron-rich due to the 
beta-captures that are favored by the increase of pressure of the degenerate elec- 


(Haensel et al. 2007. Chapt. 7). In this 
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txons. Neutrons start to drip out of nuclei at density Pdup = 4.3 x lO" g cm^^. Thus 
at p > Pdrip neutron-rich nuclei are embedded in the sea of quasi-free neutrons. 

At the bottom of the crust, the nuclei may take r odlike and platelike shapes t o 
compose so called pasta phases of nuclear matter (IPethick and Ravenhalll 1995). 
Then they form a mantle with anisotropic kinetic properties dPethick and Potekhin 


19981) . Thermodynamic stability of the pasta phase state and, therefore, the exis¬ 


tence of the mantle depends on the model of nuclear interactions. iLorenz et al.l (119931) 
de monstrated stability of the pasta phas es at p > g cm^^ for the FPS EoS model 
of Pandharipande and Ravenhall j 19891). but they were not f ound to be stable in m od- 
ern EoS models SLy ( Douchin and HaensefcOOl ) and BSk ( Pearson et al.ll2012 ). 

The strong gravity drives the rapid separation of chemical elements in the crust 
and the ocean. Estimates of characteristic sedimentation time range from s econds to 
mont hs, depending on local conditions and composition (see, e.g., Eq. 20 in IPotekhin 


20141) . which is a very short timescale compared to the stellar age. Therefore the 


envelopes are thought to be made of chemic ally pure layers, which are separate d by 
narrow transition bands of diffusive mixing ( De Blasi32000t Chang et~51l2010h . 


2.2 Thermal evolution equations 


The multidimensional hea t transport and thermal evolution equations in a locally fla t 
reference frame read (e.g., Aguilera et al.ll2008 : Pons et a02009l: Vigano et al. 11201 3 ) 


dT 


CvC^P — + V • (e'^^F) = e^^(// -Qy), F = 


(4) 


where F is the heat flux density, H is the heating power per unit volume, Cy 
is specific heat (Sects. 12.5113.21 and I4.2l i. Qy is neutrino emissivity (Sects. 12.61 
[331 US, K is the thermal conductivity tensor (Sects. 12.7113.41 and I4.41 i. and V = 
(^sin0)^^(9/(9^) in compliance with Eq. (|T]i. The inner 
boundary condition to the system of equations (|4|i is F = 0 at r = 0. The outer 
boundary condition is determined by the properties of a heat-blanketing envelope, 
which serves as a mediator of the internal heat into the outgoing thermal radiation. It 
will be considered in Sect. |5] Solutions to the thermal evolution equations and their 
implications are briefly reviewed in Sect.j^ 

Eor weak magnetic fields, we can assume that the temperature gradients are es¬ 
sentially radial, and that in most of the star volume (inner crust and core) the con¬ 
ductivity tensor is simply a scalar quantity times the identity matrix. In this limit, 
corrections for deviations from the ID approximation have little effect on the total 
luminosity. However, for strong fields and neutron stars with locally intense inter¬ 
nal heating sources, such as magnetars, a more accurate description, beyond the ID 
approximation, must be considered. 


2D calculations of thermal structur e and evolution of strongly magnetized neutron 


stars have been done bv several groups (Gennert et al.l2004, 2006; Perez-Azorin et al. 

2006; 

Aguilera et al. 20081 

Kaminker et al. 2012, 2014). In some of these works 

(Gennert et al. 

2006; 

Perez-Azorin et al. 

209^; 

Affuilera et al. 2008|). neutron-star 


models with superstrong (B^10'^-10*®G) toroidal magnetic fields in the crust were 
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considered, in addition to the less strong {B ^ 10*^ - lO'^ G) poloidal component that 
penetrates from the crust into the magnetosphere. The latter models help to explain 
the strongly non-uniform distribution of the effective temperature o ver the neutron- 
star surface and the possib le energy sourc e for magnetars ou tbursts ( Pons and Pernal 
2011; Pons and Reall2012|) . Only recently ( Vigano et al.l2013 ). the fully coupled evo¬ 


lution of temperature and magnetic held has been studied with detailed numerical 
simulations, which allow one to follow the long-term evolution of magnetars and 
their connection with other neutron star classes. Some results of such calculations 
will be discussed in Sect.|6] 


2.3 Basic plasma parameters 


In this section we introduce several basic parameters of Coulomb plasmas that are 
used below. To be concrete, we start with electrons and ions (including bare atomic 
nuclei). When other charged particles are present, their respective parameters are 
dehned analogously, with the obvious replacements of particle mass, charge, number 
density, etc. 

Since the major constituents of the neutron-star matter are mostly degenerate, an 
important parameter is the Fermi energy, which (without the rest energy) equals 


ep = c 


\J {mcY + {pvY ~ 


(5) 


where m is the particle mass, and pp is the Fermi momentum. For instance, for the 
Fermi gas in the absence of a quantizing magnetic field, p^i = fi where n 

is the number density, and h is the reduced Planck constant. It is convenient to use 
the dimensionless density parameter related to the Fermi momentum of electrons, 
Xt = ppe/i'fieC, where nie is the electron mass. In the outer core and the envelopes, as 
long as the baryons are non-relativistic, ~ (pgTg)^/^, where T, is the number of 
electrons per baryon and pg = p/10® g cm^^. 

Thermal de Broglie wavelengths of free ions and electrons are usually defined 


as Ai = InT?/m{r and Xe = \JInT?/nieT, where nii = Ania is the ion mass, and 
mu is the unified atomic mass unit. Here and hereafter, we use T in energy units 
and suppress the Boltzmann constant (i.e., 10® K = 86.17 eV). The quantum effects 
on ion motion are important either at Ai > ai or at T <C 7/, where Tj, = hcOp is the 


1 /O 

ion plasma temperature, and (Op = (A-ne^riiZ^/irii) ' is the ion plasma frequency. 
Debye temperature of a crystal ©d is closely related to the plasma temperatu re. In 
the harmonic approximation for the Coulomb crystal, ©d ~ 0.45 Tp (ICarnll96lh . 

The Coulomb plasmas are called strongly coupled if the parameter Fq = 
{Ze)^/aiT, which estimates the electrostatic to thermal energy ratio, is large. Here, 
is the ion sphere, or Wigner-Seitz cell, radius, and tii is the ion num¬ 
ber density. If the plasma only consists of electrons and non-relativistic ions of one 
kind, which is typical for neutron-star envelopes, then 


Tp = 7.832 (Z/A)v^x 10® K, Fc = 22.747Z®/3(py,)i/Vr6. (6) 
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Analogously, IJ,,, = hiAne^nelnieYI'^ = 3.34 ^P(,ZjA x 10* K is the electron 
plasma temperature. Other plasma parameters, which become important in a strong 
magnetic field, will be considered in Sect. 14.11 


2.4 Equation of state 

2.4.1 Equation of state for the outer crust and the ocean 


The composition of the outer crust and the ocean of a neutron star is particularly 
simple; their basic constituents are electrons and nuclei, which, to a good accuracy, 
can be treated as pointlik e. The EoS of such electron-ion plasmas is well known (see, 
e.g., Haensel et al.ll2007 . Chapt. 2, and references therein). 

The ions thermodynamic state will go from liquid to solid as the star cools, and in 
the solid state from a classical to a quantum crystal. It is generally assumed that the 
ions form a crystalline solid and not an amorph ous one. This assum ption is confirmed 
by molecular dynamics numerical simulations ( Hughto et al.ll2Ql 1 1 and corroborated 
by the analysis of observations of neutron-star crust cooling after an accretion episode 
(see Sect. iji. 

The simplest model of the electron-ion plasmas is the one component plasma 
(OCP) model, which considers Coulomb interactions of identical pointlike ions and 
replaces the degenerate electron gas by a static uniform charge-compensating back- 
ground. The OCP h as been studied analytically and numerically in many papers (see 
Haensel et al. 20071 Chapt. 2, for references). In the classical regime (T ^ Tp) its 
thermodynamic functions depend on a single parameter Fq. At Fq 1 the ions form 
a Debye-Hiickel gas, with increasing Fz the gas gradually becomes a liquid, and with 
further increase of Fq the OCP liquid freezes. An analysis of Monte Carlo simu- 
lations of the OCP shows tha t its ground state is crystalline when Fq > F^ — 175 
( Potekhin and Chabriei 2000h . However, supercooling cannot be excluded up to a 
value Fc cz 250. In deed, Monte Carlo simulations of freezing of classical OCP 
( DeWitt et alJ 1993 ) indicate that, as a rule, the ions do not freeze at the equilibrium 
melting temperature but form a supercooled fluid and freeze at lower T (depending 
on initial conditions and other parameters). This happens because the phase transition 
is really tiny. 

At r < Tp, the quantum effects on ion motion become significant. Then thermo¬ 
dynamic functions depend not only on Fq, but also on a. The quantum effects are 
especially important for the solid neutron star crust at high densities, although they 
can also be significant in the deep layers of the ocean composed of light elements (for 
instance, they prevent solidification of H and He plasmas). The free energy per unit 
volume of an OCP crystal can be written as 


^lat ^ I 3 ^ Eth ^ Eah 

— — fi\(~yC) -— yi\ ii] Ir) H--, 

V 2 ^ V V ' 


(7) 


where the first term is the classical static-lattice energy, Cq ~ —0.9 is the Madelung 
constant, and the next two ter ms describe t hermodynamics of the phonon gas in the 
harmonic approximation (e.g.. lKittellll963h : the second term accounts for zero-point 
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quantum vibrations, andFth/V = 'in\T (ln[l — exp(—7ia);ta/r)])ph is the thermal con¬ 
tribution. Here, u\ = {(i>ka)ph/(0^ ~ 0.5 is the reduced first moment of phonon fre¬ 
quencies cuta, and (...)ph denotes the averaging over phonon polarizations a and 
wave vectors k in the first Brillouin zone. The last term in Eq. 0 arises from an- 
harmonic corr ections, which have only been s tudied in detail in the classical regime 


{T ^ Tj,; e.g.. iFarouki and Hamaguchilll993l and references therein). An an alytical 


extrapolation of for any T was proposed in Potekhin and Chabrieil ( 20ld) . 

For mixtures of various ion species, the simplest evaluation of the thermodynamic 
functions is given by the average of their values for pure subst ances, weighed wit h 
their number fractions, which is called the linear mixing rule (iHansen et al.lll977h . 
The linear mixing rule is accurate within a few percent, if the electrons are strongly 
degenerate and fc > 1 for each of the ion species in the mixture. However, this ac¬ 
curacy may be insufficient for such subtle phenomena as melting/freezing or phase 
separatio n in the Coulomb plasmas. Corrections to the linear mixing rule were ob¬ 
tained by Potekhin et al. ( 2009 ). Medin and Curnmh^ ( 2010l) used these results to 
construct a semianalytical model for prediction of the composition and phase state of 
multicomponent mixtures. IHughto et alJ (l2012h confirmed the qualitative validity of 
this model by molecular dynamics simulations. 

The OCP model is a reasonable first approximation, but in reality the electrons 
do not form a uniform background; they interact with each other and with the ions, 
which gives rise to exchange-correlation and polarization corrections, respectively. 
The polarization corrections are appreciable even for strongly degenerate plasmas. 
For instance, they can substan tially shift the melting transition away from ic = 175 
dPotekhin and Chabrietil2013h . In the outer envelopes of a neutron star, the electron 
degeneracy weakens, and one should take the T -dependence of their EoS into ac¬ 
count. Analytical fits for all ab ove-mentioned contribu t ions to the E oS of electron- 
ion plasmas were presented by Potekhin and Chabrieil ( 2010l 2013 ). Their Fortran 
implementation is publicly available at http: //www. iof f e. ru/astro/EIP/. 

An essential input for calculating the EoS is the chemical composition of the 
plasma. The ground state of the matter in the outer crust can be found following 
' Bavm et al. ( 197 ih . The procedure, based on the minimi z ation of the 


the method of _ 

Gibbs free energy per nucleon, is described in detail in Haensel et al. ( 2007h . The 
structure of the crust is completely determined by the experimental nuclear data up 
to a density of the order p~6xl0**’g cm^^. At higher densities the nuclei are so 
neutron rich that they have not yet been experimentally studied, and the composition 
of these dense layers is model dependent. However, this model de pendence is not very 
significant in the models based on modern nu clear physics data (IHaensel and Pichon 


1994 : Riister et ^12006 : Pearson et al.ll2011|) . 


While a newly born neutron star is made of hot matter in nuclear equilibrium, 
its subsequent evolution can lead to the formation of regions where the matter is out 
of nuclear equilibrium. This is the case of an old accreting neutron star. Burning of 
the helium layer near the surface is followed by electron captures and beta decays 
in deeper layers. The growing layer of the processed accreted matter pushes down 
and eventually replaces the original catalyzed (ground-state) crust. T hese processes 
were studied by several authors (see Haensel and Zduni5ll99d 20081 and references 
therein). 
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2.4.2 Equation of state for the inner crust and the core 

The pressure in the inner crust of a neutron star is largely created by degenerate neu¬ 
trons. However, the electrons and nuclei may give an important contribution to the 
heat capacity (see Sect. 12.51) . In the core, there are contributions from neutrons, pro¬ 
tons, electrons and muons (and other particles in the exotic models, which we do not 
consider here). Different theoretical EoSs of the neutron fluid and npe^l matter have 
been proposed, based on different methods of theoretical physics: the Brueckner- 
Bethe-Goldstone theory, the Green’s function method, variation al methods, the rela- 
tivistic mean held theory, and the de nsity functional met hod (see Haensel et al .120071 


Chapt. 5, for review). The model of lAkmal et al.l (119981) (APR) has been often cited 


as the most advanced one for the core. It was derived from the variational principle of 
quantum mechanics, under which an energy minimum for the trial wave function was 
sought. The trial function was constructed by applying the linear combination of op¬ 
erators describing admissible symmetry transformations in the coordinate, spin, and 
isospin spaces to the Slater determinant consisting of wave functions for free nucle¬ 
ons. The APR EoS exists in several variants, which differ in the effective potentials of 
nucleon-nucleon interaction and in relativistic boost corrections. The potentials bor¬ 
rowed from earlier publications were optimized so as to most accurately reproduce 
the results of nuclear physics experiments. 

Many theoretical neutron-star EoSs in the literature consist of crust and core 
segments obtained using different physical models. The crust-core interface there 
has no physical meaning, and both segments are joined using an ad hoc matching 
procedure. This generally leads to thermodynamic inconsistencies. The EoS mod¬ 
els that avoid this problem by describing the core and the crust in frames of the 


(Pandharinande and Ravenhalllll989 

Lorenz et al.lll993). SLv ( 

Douchin and Haensel 

2001 

), and BSk ( 

Gorielv et alJl2010 

Pearson et alJl201 ll 2012 

) EoS families. All of 


them are based on effective Skyrme-like energy density functionals. In particular, 
the BSk21 model is based on a generalized Skyrme functional that most success¬ 
fully satisfies various experimental restrictions along wi th a number of astrophysical 
requirements (see the discussion in Potekhin et al.ll2013l) . 


2.5 Specific heat 

2.5.1 Specific heat of electron-ion plasmas 


The two components that largely dominate the specific heat Cv in the crust are the 
electron gas and the ions. In the neutron-star crust and core, the electrons form an 
ultra-relativistic highly degenerate Eermi gas, and their contribution in the heat ca¬ 
pacity per unit volume is simply given by 

^2 


pie 

3h^c 


T = np 


-T. 


PVeC 


( 8 ) 


In the ocean, where the density is lower, approximation (|3) may not work. Then it is 
advisable to use accurate approximations, cited in Sect. 12.4.11 
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Fig. 1 Left panel: Heat capacity per ion versus T (bottom axis) and Ic (top axis) for ®^Ni at p — 10*^ 
g cm“^. The solid line displays the total normalized heat capacity Cyjnw the long-dashed line shows this 
quantity for a classical Coulomb lattice of ions, including harmonic and anharmonic terms; the short- 
dashed line is the harmonic-lattice approximation in the solid phase; the dot-dashed line is the same plus 
anharmonic and electron polaiization corrections in the solid phase. The dotted line is the electron Fermi 
gas contribution. The vertical line is the OCP melting point Ic = 175, and the hatched band shows the 
range Ic = 150-200, where melting is expected to occur in realistic conditions. The inset illustrates the 
competition between the electron and ion contributions t low T. Right panel: Heat capacity per baryon as 
function of mass density from the ocean throughout the crust and core of a neutron star at T = 10^ K. 
The solid line displays the total Cv/^b’ ^nd the other lines show its constituents due to the electrons {e~), 
neutrons (n) in the inner crust and core, nuclei (A^), including electrostatic terms in the ocean and crust but 
neglecting the neutron entrainment effects in the inner crust (Sect. 1232), protons {p) and muons (p ) in 
the core, assuming that the nucleons are non-superfluid. For comparison, the thick long dashes display an 
example of the total Cv/^b in the inner crust and core in the case of superfluid nucleons (Sect. [22). The 
top axis shows the volume contained inside a sphere with given p for a 1.4 Mq neutron star. The stellar 
structure and composition con'espond to the BSk21 EoS model. 


In Fig. [T] we show the temperature and density dependences of the normalized 
heat capacity of the ground-state matter in a neutron star. The left panel illustrates 
the dependence of Cv/n^ on T, and the right panel the dependence of Cy/nb on p. 
Since the electron polarization effects shift the melting temperature (Sect. imi l. the 
phase transition may occur anywhere within the hatched region around the vertical 
line fc = 175 in the left panel. 

When the temperature of the Coulomb liquid decreases, the heat capacity per ion 
increases from the ideal-gas value Cv,i/ni = | at T ^ I'm to, approximately, the sim¬ 
ple harmonic lattice value Cv,i/ni = 3 at T < Tn, (the Dulong-Petit law for a classical 
harmonic crystal). This gradual increase is due to the Coulomb non-ideality in the 
liquid phase, which effectively smears a p hase transition bet ween the strongly cou¬ 
pled Coulomb liquid and OCP crystal (see Baiko et ^ 1998h . With further cooling, 
quantum effects suppress the heat capacity. Once the crystal is deep into the quantum 
regime its specihc heat is given by the Debye result 


(D) 

<7 = Hi 


127t:^ 


5 


T 

©D 


3 


( 9 ) 
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The calculations of Baiko et al. ( 2001bl) show that the Dulong-Petit law applies at 
temperatures down to T ~ Tp, while the Debye value of Eq. (|3) is attained when 
T < 0.1 Tp. The same authors present a simple analytical approximation for the heat 
capacity of a harmonic Coulomb crystal, accurate to a few parts in 10^ at any T. 

However, the harmonic OCP model is an idealization. The anharmonic correc¬ 
tions and electron polarization cori'ections (Sect.|2ATJ can amount up to (10-20) % 
of Cy.i. Because of the anharmonic effects, Cv,i/ni is not equal to 3 exactly, but is 
^ 10% larger ?AT — Tm- If the above-mentioned supercooling takes place in stellar 
matter, various fluid elements solidify at different T below Tm, and the average heat 
capacity, as function of te mperature, can conta in a bump, associated with latent heat 
releases (see Sect. 2.4.6 of Haensel et al.lE007l for a discussion). 

The right panel of Fig. [T]shows the density dependence of the total heat capacity, 
normalized per one nucleon, Cy/ub, throughout the neutron star from the ocean to 
the core, and partial co ntributions to Cy/ny,- Different particle fractions are adopted 
from the BSk21 model ( Gorielv et al.ll201ot Pearson et al. 201 li 2012 ). as fitted by 
Potekhin et al.l (120 13h . Here, we have mostly neglected the effects of nucleon super¬ 
fluidity to be discussed in Sect. [3] The importance of these effects is demonstrated, 
however, by the heavy long-dashed line, which displays the total normalized heat 
capacity suppressed by nucleon superfluidity (see Sect. l3.21 l. 


2.5.2 Specific heat of neutrons 


In the inner crust, besides electrons and nuclei, there are also neutrons. In a thin layer 
at densities p just above the neutron drip point Pdrip^ the dripped neutrons are not 
paired (non-superfluid) and largely dominate Cy. Heat capacity of strongly degener¬ 
ate non-superfluid neutrons can be accurately evaluated using the above-referenced 
analytical fits, but since the neutrons are strongly degenerate almost everywhere in 
the neutron star, the simpler Sommerfeld result for Fermi gases at T -C Ep is usually 
applicable, 

Cy,x ^ ~Z ) ( 10 ) 

Ep.x 

where x stands for the fermion type (x = n, p, e, p). For neutrons at p only slightly 
above PdUp, however, the l atter formula is ina ccurate because gp „ is not suffi¬ 


ciently large. For this reason, Pastore et alJ (l2015h proposed an interpolation between 

iu, 


Eq. (ITOl) and the ideal-gas limit Cy = 




( 11 ) 


They also showed that co rrections due to the coupling to phonons (e.g.. Sect. 1.4.4 in 
Bavm and Pethicklll99lh turn out to be unimportant for Cy^„. Approximation (fTTT) is 
accurate within 17% for non-relativistic Fermi gases at any density. For a relativistic 
Fermi gas, we can preserve this accuracy by using Eq. (|5]) for Ep x and multiplying 
both Tci and prefactor | by the ratio {mc^ -f lOT)/(mc^ + 5T). 

With further density increase, the neutrons become superfluid (Sect. |3, and then 
their contribution to Cy nearly vanishes. H owever, even in a superfl uid state, the neu¬ 
trons have a dramatic effect on Cy. Indeed. Flowers and Itohl ( 1976 ) noticed that since 
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Fig. 2 Iso-contour lines of Cy in the crust, labeled by the value of log£Q(cv/ergcm ^ K ^). Also shown 
are the melting curve and the critical temperature for neutron ^So superfluidity, Tq, the electron and 
ion plasma temperatures, and 7p respectively, the Debye temperature, 6 )d — 0.45rp, that marks the 
transition from classical to quantum solid and O-Up be low which the wholly quan tum crystal regime is 
realized. The outer crust che mical composition is from iHaensel and PichonI j 19941) and inner crust from 
iNegele and Vautherinl jl973l) with the neutron drip point at p^j-jp — 4.3 x 10^* g cm The electron con¬ 
tribution dominates in the two dark-shadowed (green) regions at high T and p below Pdrip and at low T 
and high p, while neutrons dominate in the light-shadowed (yellow) region at high T and p above Pdrip. 
and ions dominate in the intermediate regime. The right panel only displays the inner crust but assuming 
that about 80% of the dripped neutrons are entrained, illustrating the resulting increase in Cy i, mainly due 
to the strong reduction of Tj, and 0 d, significantly extending the regime where Cyp dominates over Cy^^,. 


free neutrons move in a periodic potential created by lattice of atomic nuclei, their 
energy spectrum should have a band structure, w hich can affect kinetic and neutrino 
emission phenomena involving the free neutrons . IChamell (l2QQ5h calculated the band 
structure of these neutrons, in much the same way as electron band structure is calcu¬ 
lated in solid state physics. The effect of this band structure is that a large fraction of 
the dripped neutrons are “locked” to the nuclei, i.e., the thermal motion of the nuclei 
entrains a significant part of the dripped neutrons resulting in a strongly increased 
ion effective mass mgffj. This increase wieir,i significantly increases Cyj in the 


quantum regime since c 


(D) 


( Charnel et alj|2013l) . 


The overall “landscape” of crustal specific heat is illustrated in Fig.|2l For highly 
degenerate electrons Cv,e T, while for ions Cv,i decreases as according to Eq. (|3), 
therefore the electron contribution dominates at T <C ©d. and the ion contribution 
prevails at T > ©d (cf. the inset in the left panel of Fig. [T]|. On the other hand, in the 
non-degenerate regime Cv,e/cv4 ~ Z, therefore the contribution of the electrons dom¬ 
inates again for Z > 1 at T > ep in the liquid phase (also cf. the left panel of Fig. [T]i. 
The effect of dripped neutron band structure on low-level collective excitations in the 
inner crust and the resulting increase of Cv,i is illustrated in the right panel of Fig. |2] 
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2.5.3 Specific heat of the core 

The specific heat is simpler to evaluate in the core than in the crust but it has larger 
uncertainties. The core is a homogeneous quantum liquid of strongly degenerate 
fermions, and its specific heat is simply taken as the sum of its components con¬ 
tribution: Cy = LxCv,x where x stands for neutrons (n), protons (p), electrons (e), 
muons (p), and any other component as hyperons or quarks that may appear at high 
densities. For each fermionic component, one can use Eq. cnii, but for baryons one 
should replace the bare fermion mass iWx by an effective mass m*, which encapsu¬ 
lates most effects of interactions. In principle, m* should be calculated from the same 
microphysical interaction as employed for the EoS; cf. Sect. 12.6.31 Eor leptons (e 
and p), interactions have a negligible effect on m* and the bare fermion mass value 
can be used. The nucleon heat capacity in the core is strongly affected by pairing 
(superfluidity) effects, as discussed in Sect. 13.21 


2.6 Neutrino emissivity 


The neutrino luminosity of a neutron star is, in most cases, strongly dominated by 
the core contribution, simply because the core comprises a lion’s share of the total 
mass. The crust contribution can, however, prevail in the case of strong superfluidity 
in the core, which suppresses the neutrino emissivities. Crust neutrino emission is 
also essential during the early thermal relaxation phase (the first few decades of the 
life of the star), or the crust relaxation after energetic transient events (e.g., strong 
bur sts of accretion i n X-ra y binaries and flares in magnetars). 


Yakovlev et al.l (1200 Ih reviewed the main neutrino emission mechanisms in neu¬ 


tron star crusts and cores and collected fitting formulae for the neutrino emissivity 
in each reaction as a function of density and temperature. The summary of the most 
important processes is given in Table [1] The last column of this table contains refer¬ 
ences to the analytical fitting formulae that can be directly employed to calculate the 
relevant emission rates. These processes are briefly described below. 


2.6. J Neutrino emission in the crust 

There is a variety of neutrino processes acting in the crust. In a non-magnetized crust 
the most important ones are the plasmon decay process and the electron-ion brems- 
strahlung process (see Table[T]l. The pair annihilation process can be also important 
if the crust is sufficiently hot. 

The total emissivity from the sum of these processes is illustrated in the left panel 
of Eig. [3 The first thing to notice is the enormous range of values of Qv covered in 
the p — T range displayed in this figure, which spans 26 orders of magnitude. This is 
a direct consequence of the strong T dependence of the neutrino processes. The pair 
annihilation process is efficient only at low densities and very high temperatures, but 
when T ^ Ip e very few positrons are present and the process is strongly suppressed. 
In the whole range of this plot, rF,e ^ 10^*^ K but pair annihilation still dominates 
at low p and high T. In the opposite high-p and low-T regime the dominant process 
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Table 1 Main neutrino emission processes^ 


Process / Control function 

Symbolic notation^ 

Formulae for Qy and/or R 




In the crust 


1 

Plasmon decay 

r ^v + v 

Eqs. (15)-(32)of [1] 

2 

Electron-nucleus bremsstrahlung 

e' 

^+N^e^+N+v + v 

Eqs. ( 6 ), (16)-(21)of [2] 

3 

Electron-positron annihilation 

e' 

“ 4- e+ ^ V + V 

Eq. (22) of [3] 

4 c 

Electron synchrotron 

e' 

^ Ae“+v + v 

Eq. (48)-(57) of [3] 




In the core 


!■* 

Direct Urea 
(Durca) 

pFe + Ve , 
p-\-e~ — )• « + Vg 

Eq. (120) of [3] 


Magnetic modification^ 

R 

(D) 

6 

Eqs. (247)-(250) of [3] 


Reduction factors^ 

R 

(D) 

X 

Eqs. (199), (202)-(206) of [3] 

2 

Modified Urea (Murca) 

(neutron branch) 

n-\-n ^ n + p + e +Ve, 
n-\- p-{-e~ — >n + n + Ve 

Eq. (140) of [3] 


Reduction factors^ 

R 

(Mn) 

X 

Appendix of [41 

3 

Murca 

p + n^ p + p + e +Vc, 

Eq. (142) of [3], con'ected at 


(proton branch) 

p + p + e^ p + n + Ve 

3pFp > PVn + PFe as pet [4] 


Reduction factors^ 

R 

(Mp) 

X 

Appendix (and Eq. (25)) of [41 

4 

Baryon-baryon 

bremsstrahlung 

\ 

rt + n— >rt + n + v + v 
« + ^H+/? + V + V 

P+P^ p+p+v+v 

Eq. (165) of [3] 

Eq. (166) of [3] 

Eq. (167) of [3] 





Eqs. (221), (222), (228) of [3] 




and Eq. (60) of [4] 


Reduction factors^ 

< 


Eq. (220), (229) of [3] 





and Eq. (54) of [4] 

Eq. (221) of [3] 

S" 

Cooper pairing of baryons 

\ 

' n + n^ [nri\ + v + v 
+ M + V + V 

Eqs. (236), (241) of [31, 
corrected as per [5] (Sect.fT^ 

gc.e 

Electron-fluxoid 

bremsstrahlung 

e' 

^-l-/->e“+/+v + v 

Eqs. (253), (263), (266)-(268) 
of [3] 


Notes. ^ References: [11 [Kanto r and Gusakov! j 2007l): [2] lOfengeim et alj j2Q14l) : [3] lYakovlev et al.l 
j 20 Qlh : [4] iGusakovI j 20 Q 2 h : [51 iLeinsonI j 2 QQ^l 2 QlQh . ^ F means a plasmon, e~ an electron, e'^ a 
positron, v a neutrino, v an antineutrino (in general, of any flavor, but Ve or stands for the electron 
neutrino or antineutrino, respectively), p a proton, n a neutron, [pp] and [nn] their paired states, N stands 
for an atomic nucleus, and / for a proton fluxoid. At densities where muons are present, they participate in 
the Urea and bremsstrahlung processes fully analogous to the processes 1, 2, 3, 6 in the core (see details in 
Ref. [11). R with subscripts/superscripts signifies a control function (correction factor) due to supeidluidity 
or magnetic field. Subscript x in Rx substitutes for different superfluidity types (proton or neutron, singlet 
or triplet); B indicates magnetic field. ^ The effect of strong magnetic field (see Sect. IS). ^ At densities 
beyond the Durca threshold (see Sect. 1122). >= The effect of superfluidity (see Sect. IS). 


is electxon-ion bremsstrahlung, for which oc At intermediate T and p 

the plasmon decay process is most important and, when it strongly dominates, its 
emissivity behaves as qF 

The right panel of Fig.[3]illustrates the density dependence of 
in either ground-state or accreted crust of a neutron star with T = 10^ K. Pair an- 
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Fig. 3 Neutrino emissivity Qy in a non-magnetized crust from the pair annihilation, plasmon decay, 
and electron-ion bremsstrahlung processes. Left panel: The contour lines are labeled by the value of 
logio[2v/(ei'gcrn^^s“*)]. Regions where the pair, plasma, and bremsstrahlung processes dominate are 
indicated: the boundaries happen to be quite well described by the two dotted lines that show j Tp,, and 
j^Tp e. (Also indicated is the ion melting curve, dashed line.) Right panel: Density dependences of Qy for 
the ground-state nuclear matter (solid lines) and for the accreted crust (dashed lines) at T = 10^ K. The 
dotted line represents an older fit to the bremsstrahlung process (see text for detail). 


-)(pl) 


nihilation is negligible in this case. Qy ^ is calculated according to lYakovlev et al 
( 2001 ) and according to Ofengeim et al.l ( 20141) . For comparison, an older 

fit to ( Kaminker et al 


1999h is plotted by the dotted line. Th e ground-state 


composition and the nucle ar size are d escribed by the BSk2 1 model (iGorielv et al 


2010t Pearson et al. [201 ^ as fitted by Potekhin et al. 2013). The accreted compo- 
sition is taken fromIHaensel and Zdunik ( 199ol) : in~this case the approximation by 
Itoh and Kohvamal (Il983h is used for the nuclear size. 


The band structure of the energy spectrum of neutrons in the inner crust, which 
was mentioned in Sect. 12.5.11 should reduce the neutrino reactions of the brems¬ 
strahlung type and initiate an additional neutrino emission due to direct inter-band 
transitions of the neutrons, in analogy with Coop er pairing of neutrons discussed in 
Sect. 13.31 These effects have been mentioned by lYakovlev et alJ (l20Qlh . but remain 
unexplored. 

Electron and positron captures and decays by atomic nuclei (beta processes), 
which accompany cooling of matter and non-equilibrium nuclear reactions, pro¬ 
duce neutrino emission. A pair of consecutive beta capture and decay reactions is 
a nuclear Urea process. Urea processes involving electrons were put forward by 
Gamow and Schoenberg! (Il94lh . while those involving positrons were introduced by 
PinaevI (119641) . In the neutron star crust, the appropriate neutrino luminosity de¬ 
pends on cooling rate and should be especially strong at T ^(2^) xlO^ K when 
the main fraction of free neutrons is captured by nuclei. However, there are other ef¬ 
ficient neutrino reactions open at such te mperatures, which ma ke the neutrino emis¬ 
sion due to beta processes insignificant ( Yakovlev et al. 200 ih . On the other hand, 
heating produced by non-equilibrium nuclear reactions (the deep crustal heating. 
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Haensel and Zdunik 1990l) that accompany accretion mentioned in Sect. 12,4.11 


may 

be more important than non-equilibrium neutrino cooling. _ 

There are a number of other neutrino-emission processes ( Yakovlev et alJ 2001 ). 
which are less efficient than those listed in Table [T] In the inner crust with dripped 
neutrons, n — n bremsstrahlung is very efficient but it is suppressed by pairing and, 
hence, only acts in the layers where T > Tc, where Ten is the neutron pairing crit¬ 
ical temperature (Sect. [IB- This process operates in a wide range of densities and 
temperatures, and the density dependence of its emissivity is generally smooth. Neu¬ 
trino emission from the formation and breaking of Cooper pairs makes a significant 
contribution, much stronger than the bremsstrahlung, but is confined to a restricted 
region of p and T (Sect. ISJl l. In the presence of a very strong magnetic field, some of 
the above-mentioned processes are modified, and new channels for neutrino emission 
may open (Sect. nil. 


2.6.2 Neutrino emission in the core 


Yakovlev et al.l ( ^Qlh discussed a wealth of neutrino reactions which may be impor¬ 
tant in the npejiAL^ matter in a neutron star core, which include 


1. 8 direct Urea (Durca) processes of the electron or muon production and capture 
by baryons (baryon direct Urea processes), 

2. 32 modified Urea (Murca) processes, also associated with the electron or muon 
production and capture by baryons (baryon Murca processes), 

3. 12 processes of neutrino-pair emission in strong baryon-baryon collisions (baryon 
bremsstrahlung), 

4. 4 Murca processes associated with muon decay and production by electrons (lep¬ 
ton Murca process), 

5. 7 processes of neutrino pair emission in Coulomb collisions (Coulomb brems¬ 
strahlung). 


In this paper we basically restric t ourse lves to the npep matter. We refer the reader 
to the review by [Yakovlev et al.l OOOll) for the more general case, as well as for a 
discussion of other exotic models (such as the pion or kaon condensates). It appears 
that the reactions that proceed in the npep matter are often sufficient for the neutron- 
star cooling, even when the appearance of the and A hyperons is allowed. The 
reason is that these hyperons can appear at high densities only, where the baryon 
Durca processes are likely to be allowed and dominate, for realistic EoSs. 

The Durca cycle consists of the beta decay and electron capture processes (see 
Table[B- They are threshold reactions open at sufficiently high densities, and not for 
every EoS model. Eor the degenerate nucleons they are only possible if the proton 
fraction exceeds a certain threshold. In the npe matter (without muons) this threshold 
is Ri 11%, which follows readily from the energy and momentum conservation com¬ 
bined with the condition of electric charge neutrality of matter. Indeed, for strongly 
degenerate fermions the Pauli blocking implies that the reaction is possible only if 
the energies of the reacting particles are close to their respective Eermi energies. 
Then the momentum conservation assumes the inequality ppn ^ PFe + PFp, that is 
tin ^tie +Hp . Eor the npe matter Mg = because of the charge neutrality, there- 





















18 


A.Y. Potekhin, J.A. Pons, D. Page 


fore n„ ^ 8np, or Hp ^ nb/9, where nb is the total baryon number density. The pres¬ 
ence of muons can increase this threshold by several percent. If ^ ppn — ppp, then 
the muon Durca process adds to the electron Durca. 

If allowed, the Durca processes produce a rapid (enhanced) cooling of neutron 
stars. If they are forbidden, the main reactions are those of the baryon Murca and 
bremsstrahlung processes which produce a slow (standard) cooling. The Murca pro¬ 
cess is a second order process, in which a bystander neutron or proton participates 
to allow momentum conservation (see Table [TJ. Since this process involves five de¬ 
generate fermions, instead of three for the Durca process, its efficiency is reduced, 
simply by phase space limitation, by a factor of order (T/ep)^, which gives an over¬ 
all temperature-dependence instead of T®. This reduction, for typical conditions 
in the neutron-star core, amounts to 6 orders of magnitude. It is certainly the domi¬ 
nant process for not too high densities in absence of pairing, and is the essence of the 
“standard cooling scenario”. However, in presence of superfluidity, neutrino emission 
by the formation of Cooper pairs (Sect. 13.31 1 can dominate over the Murca process. 

Other neutrino reactions in the core involve neutrino-pair bremsstrahlung in Cou¬ 
lomb collisions lepton modified Urea processes, electron-positron annihilation, etc. 
All of them are not significant under the typical conditions in the non-exotic core. 
For instance, the plasmon decay process that is efficient in the neutron star crust 
(Sect. 12.6. Il l is exponentially suppressed in the core, because the electron plasmon 
energy in the core (^ /icUp ^ 10 MeV) is much larger than the thermal energy. 

In a strong magnetic field penetrating into the core, some of the above-mentioned 
processes can be modified, and new channels for neutrino emission may open (see 
Sect.lH). 


2.6.3 Remarks on in-medium effects 


Neutrino emissivity Qy may be strongly modified by in-medi um (collective) effects 
at the high densities of neutron stars (see Voskresenskvll20Ql . for a review). For in¬ 
stance, these effects may result in renormalization of electroweak interaction param¬ 
eters. Moreoyer^JheJn-mediumdfectsmay open new channels for neutrino emis¬ 
sion. Voskresensky and Senatorov! (Il986h found that the direct and modified Urea 
processes appreciably exceed the estimates obtained neglecting the collective ef¬ 
fects, provided the density is sufficiently large. On the other hand, the in-medium 
effects suppress the n n bremsstrahlung in the neut ron-star core by a fac tor of 10-20 
( Blaschke et al. 1995h . According to the study by Schaab et al. (1997), the medium 
effects on the emissivity of the Murca process cause a more rapid cooling than ob¬ 
tained for the standard scenario and result in a strong density dependence, which 
gives a smooth crossover from the standard to the enhanced cooling scenario (see 
Sect. 16.Ill for increasing star masses. 

The problem of calculation of the in-medium effects in the neutron star matter is 
complicated. Various theoretical approaches were used to solve it, res ults of different 


techn iques being different typically by a factor of a few (see, e.g.. iBlaschke et al 


19951 and references therein). The renormalization of the electroweak coupling is 
usually taken into account in an approximate manner by replacing the bare baryon 
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masses mg with effective ones, m*g (e.g.. lYakovlev et al.ll2001h . The values of these 
effective masses should be taken from microscopic theories. 

The in-medium effects are also important for the Cooper-pairing neutrino emis¬ 
sion mechanism related to baryon superfluidity, as discussed in Sect. [33] below, for 
heat capacity (Sect. lZ33]) . and for baryon heat conduction (Sect. [2331) in the core of 
a neutron star. 


2.7 Thermal conductivity 

The most important heat carriers in the crust and ocean of the star are the electrons. 
In the atmosphere, the heat is carried mainly by photons. In general, the two mech¬ 
anisms work in parallel, hence K = (Q + fQ, where fQ and (Q denote the radiative (r) 
and electron (e) components of the thermal conductivity K. The radiative transfer is 
considered in Paper I. In this paper we will pay most attention to the electron heat 
conduction mechanism. Both the electron and photon heat conduction are affected by 
strong magnetic fields. We will consider these effects in Sect.|4| 

The elementary theory in which the effective collision rate v of heat carriers with 
effective mas s m* and number density n does not depend on their velocity, gives 
( Ziman|[l960h 

K — anT/m*V, (12) 

where a is a numerical coefficient; a = 3/2 for a non-degenerate gas, and a — n'^/3 for 
strongly degenerate particles. (We remind that we use energy units for T ; otherwise 
a should be multiplied by the squared Boltzmann constant.) 

The most important heat carriers and respective scattering processes that control 
the thermal conductivity K are listed in Table |2| and briefly discussed below. The 
last column of the table contains references to either analytical fitting formulae or 
publicly available computer codes for the evaluation of K. Figure 01 illustrates the 
magnitude of (Q and characteristic temperatures in the crust. 


2.7.1 Heat conduction in the outer envelopes 


Electron heat conduction is the most important process in the neutron star envelopes 
that determines thermal luminosity of neutron stars. In this case, m* = nie -^/l +xj in 
Eq. (fT2l l. and v = Vg is mostly determined by electron-ion (ei) and electron-electron 
{ee) Coulomb collisions. In the crystalline phase, the electron-ion scattering takes 
the f orm of scatter ing on phonons (collective ion excitations). The Matthiessen rule 
(e.g., Ziman 1960h assumes that effective frequencies of different collisions simply 
add up, i.e., Vg = Vgj + Vgg. This is strictly valid for extremely degenerate electrons 
( Hubbard and Lamr)ell969 ). In general case it remains a good estimate , because Vgi -f 
Vee < Vg < Vgi -f Vgg -f 5v, where 5v ^ min(Vgi, Vgg) (IZimanl 1960l) . The relative 
importance of the different types of collisions and practical formulae for evaluation 
of Vg can be different, depending on the composition and phase state of the plasma 
(see AppendixiAl. 
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Table 2 Main contributions to thermal conductivity 


Conduction type and regime 

References^ 


Photon conduction 

Eqs. (14)-(20) of[l] 


- plasma cutoff correction 

Sect. 3.3 of [2] 


- magnetic field modifications'^ 

Eqs. (21)-(23) of [1] 

2 

Electron conduction in the ocean and the crust: 

see Appendix!^ 


- Electron-ion / electron-phonon scattering 

[3] (theory), [4] (public code) 


- the effects of magnetic helds'^ 

[5] (theory), [4] (public code) 


- the effects of finite nuclear sizes in the inner crust 

[6] (theory), [4] (public code) 


- Electron scattering on impurities in the crust 

- Electron-electron scattering: 

see Appendix [A^ 


— strongly degenerate electrons 

Eqs. (10), (21)-(23)of [7] 


— arbitrary degeneracy 

see Appendix [A^ 

3 

Baryon conduction in the core 

Eqs. (7), (12), (21), (28)-(30) of [8] 


- Effects of superfluidity® 

Eqs. (45)-(48), (50)-(53) of [8] 

4 

Lepton conduction in the core 

Eqs. (4)-(6), (16), (17), (33)-(37) of [9] 


- Effects of superfluidity® 

Eqs. (45), (54)-(61), (84)-(92)f of [9] 


Notes. “References: [1] iPotekhin and Yakovle\^ i200lh: [21 iPotekfiin et alj j200^: [31 iPotelAin et alj 
(l9 99t); [4] http://www.io ffe.ru/astro/conduct/; 51 |Potekhinl jl99^~[T99 ^: [61 |Gr iedin et alJ i200lh : 
rvi lshtemiii and YakovlevI i2006h : [Sl fBaiko et al i2001ah : r9l lShtemin and Yak ovlev! l2 007h. ‘’ F or fully 
ionized atmospheres only. For partially ionized atmospheres, see references in IPotekhinl i20T3l . “ See 
Sect. 14.4.11 See Sect. 14.4.21 See Sect. Il4l f The power index 2 should be suppressed at its first 
occurrences in the third and fourth lines of Eq. (92) of Ref. [9]. 


Chugunov and Haensell (l2007h considered an alternative heat transport by the 


plasma ions (phonons in the solid OCP), which works in parallel with the trans¬ 
port by the electrons. The ion (phonon) heat conduction is usually unimportant in 
neutron stars. Although the ion thermal conductivity can be larger than the electron 
conductivity across the strong magnetic field, the multidimensional modeling shows 
that in such cases the heat is mainly transported by the electrons non-radially (i.e., 
not straight across the field lines; see Sect. IS- 


2.7.2 Heat conduction in the inner crust 


The inner crust of a neutron star is characterized by the presence of free neutrons. This 
has two important consequences. First, heat transport by neutrons can compete with 
the transport by the electrons and phonons. Second, electron-neutron scattering adds 
to the other electron scattering mechanisms considered above and in AttpendixlATl 
The thermal conductivity by neutrons, Kn, was studied in several papers (e.g.. 
Flowers and Itohlll976t [Bisnovatyi-Kogan and Romanovalll982h . A general expres¬ 
sion for Kn in non-superfluid matter is given by Eq. ([nil with n — n„, the number 
density of neutrons, m* = m*, the neutron effective mass modified by medium effects, 
and v„ = V„i -(- V„„. The neutron-neutron collision frequency, v„„, can be calculated 
in the same manner as in uniform matter of a neutron-star core (Sect. I2.7.3I I. How¬ 
ever, for strongly degenerate neutrons these collisions are much less efficient than 
the neutron-ion ones. Therefore, one can set v„ ~ V„i, at least for order-of-magnitude 
estimates. For the scattering of the neutrons by uncorrelated nuclei, v„i = ni\p„Sni, 
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p [gcm-3] 


p [gcm-5] 


Fig. 4 Iso-contour lines of the electron thermal conductivity Ke in th e crust, labeled by the value 
of logjQ(K^/ergs“^ cm“^ K“^), using the results of iGnedin et alj j20Qll) . Also shown are the melting 
curve Tm, the electron and ion plasma temperatures, Tp^e and Tp respectively, the Debye temperature, 
0D — 0A5Tp, that marks the transition from classical to quantum solid and O.lTp below which the wholly 
quantum crystal regime is realized. The crust composition is the same as in Fig. [2] The right panel only 
displays the inner crust but assuming that about 80% of the dripped neutrons are entrained: the strong 
reduction of Tp and 0 d pushes the onset of the wholly quantum regime to lower T. The dashed contour 
lines illustrate the reduction of from impurity scattering, assuming an impurity parameter Gimp = 1- 


where Vp„ = p-pn! + {PFn/tKiY '^^e neutron Fermi velocity and is the trans¬ 
port cross section. For a crude es timate at sufficiently low neutron energie s in the neu 
tron star crust one can set (e.g.. lBisnovatvi-Kogan and Romanovajll982) S„\ = 7tRl, 


where /?„ is the neutron radius of an atomic nucleus (fitted, e.g., in [Potekhin et al 


2013h . Estimated in this way, «•„ is negligible, being at least two orders of magnitude 
smaller than jq in the entire inner crust at F < 10^ K. However, v„i can be strongly 
affected by ion-ion correlations and by superfluidity (Sect. [3^ . 


In addition, the electron conduction in the inner crust is affected by the size of a 
nucleus, which becomes non-negligible compared to the mean distance between the 
nuclei, so that the approximation of pointlike scatterers is not applicable anymore. 
Then one should take into account the form factor, which depends on the size and 
shape of the charge distribution in a nucleus. A finite charge distribution reduces Vgi 
with respect to the model of a pointlike charge, thereby increasing the conductivity 
( Gnedin et al.ll2Q01 ). The effect mainly depends on the ratio of the root mean square 
charge radius of a nucleus Rch to the Wigner-Seitz cell radius ai. Gnedin et al. ( 200 lb 
presented fitting formulae for the dependences of the thermal and electrical conduc¬ 
tivities on the parameter ;iCnuc = \/5/3Rch/cii. The latter param eter has been fitted as 
function of density for modern BSk model s of nuclear matter (iPotekhin et al.ll2013h 
and for some other models (Appendix B in Haensel et alJl2007 ). 
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2.7.5 Heat conduction in the core 


The first d etailed studies of the kin etic coefficients in neutron star cores were per¬ 
formed by Flowers and Itoh ( 19791) . who constructed the exact solution of the mul¬ 
ticomponent system of transport equations in the npe matter. But since the pro¬ 
ton fraction is small and the electron-neutron interaction is weak, the kinetic coef¬ 
ficients can be split in two almost independent parts - the neutron kinetic coeffi¬ 
cients mediated by nucleon-nucleon collisions and electron kinetic coefficients me¬ 
diated by the collisions between charged particles; the proton kinetic coefficients are 
small. In the non-superfluid npepL matter, the neutrons are the main heat carriers at 
T 5,10^ K, while the heat transport by leptons a nd is competitive at T > 10^ K 
( Shternin and Yakovlevl2007t Shternin et al.ll2013 ). 


Baryon heat conduction. Flowers and Itohl ( 19791) based their calculations on the free 
nucleon scattering amplitudes, neglecting the Fermi-l iquid effects and nuc leon many- 
body effects. Their results were later reconsidered by Baiko et ^ ( 2001a ). 

The thermal conductivity is written in the form analogous to Eq. (fT^ : 


K„ = 


It ^TneXn 
3mt 


n^TnpXp 

3m* ’ 


(13) 


where the effective relaxati on times and x„ ar e provided by solution of the system 
of algebraic equations (e.g., Shternin et al.ll2013l) 


j=n.p 


Vii = 


64m* m*^T^ 


{ij = n,p), 


(14) 


where Vij are effective collision frequencies, m^ is the bare nucleon mass in vacuo, 
and Sij are the effective cross-sections. 

Many-body effects in the context of transp ort coefficients of pure neutron matter 
were first addressed bv lWambach et al.l (119931) and later reconsidered in many papers. 
There are two kinds of these effects; the three-body part of the effective potential for 
the nucleon-nucleon interactions and the in -medium effects (cf. Sect. l2.6.3l l that affect 
nucleon-nucleon scattering cross-sections. iBaiko et al.l (12001 ah calculated Sij in the 
approximation of pairwise interactions between nucleons with appropriate effective 
masses, using the Bon n potential model for the elastic nucleon-nucleon scattering 
( Machleidt et al. 1987h with and without the in-medium effects. They presented the 


results in the form Sij = where S]y corresponds to scattering of bare particles, 

and Kij describes the in-medium effects. They also constructed a simple analytical fits 
to their results for and Ka (referenced in Table|2]i. 

■ I ‘J 


r(0) , 


Shternin et al.l (12013h studied the many-body effects on the kinetic coefficients of 


nucleons in the npep matter in beta equilibrium using the Brueckner-Flartree-Fock 
(BHF) method. According to this study, the three-body forces suppress the thermal 
conductivity. This suppression is small at low densities but increases to a factor of 
^ 4 at the baryon number density of n\j = 0.6 fm^^. However, the use of the effec¬ 
tive masses partly grasps this difference. For this reason it proves to be sufficient to 
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multip ly the conductivities obtained in the effective-mass approximationjBaiko,etaJ 


200 lal) by a factor of 0.6 to reproduce the BHF thermal conductivity (IShternin et al 


2013h with an accuracy of several percent in the entire density range of interest. 


Lepton heat conduction. The up-to-date electron and muo n contributions to ther¬ 
mal co nductivities of neutron star cores were calculated by IShternin and Yakovlev 
(120071) . Their treatment included the Landau damping of electromagnetic interac- 
tions owing to the exchange of transverse plasmons. This effect was studied by 
Heiselberg and PethickI (119931) for a degenerate quark plasma, but was neglected 
in the previous studies of the lepton heat conduct ivities in the npep matter (e.g., 
Flowers and ltohiri981 ; Gnedin and Yakovlevlll995 ). 


The electron and muon thermal conductivities are additive, = Ke + and 


can be written in the familiar form of Eq. (fTSI) : 


Kp = 


7t^ Tne 


3mt 


K Tnp^lfi 

3m* ' 


(15) 


where Kg and )f^ are the partial thermal conductivities of electrons and muons, re¬ 
spectively; Mg and are number densities of these particles, m* and mjj are their 
dynamical masses at the Fermi surfaces, determined by their chemical potentials. 
In neutron star cores at beta equilibrium these chemical potentials are equal, there¬ 
fore m* = m*,. The effective electron and muon relaxation times can be written as 


( Gnedin and Yakovlevll 1 995l) 


Vu-v: 


Tg = 


e^l 


Vg-v: 


VgV^-vL,v; 


Xu = 


ne 


efi fie 


VgV^-v;„v; 


(16) 


eji fie 


yefi ' 


Vgn and Vu = I,- v„,- = v. 




yiie ■ 


V^ip are the to- 


where Vg = Vei = Vg, 
tal effective collision frequencies of electrons and muons with all charged particles 
i = e,lJ.,p', Vei and are partial effective collision frequencies, while and v^g 
are additional effective collision frequencies, which couple heat transport of the elec¬ 
trons and muons. All these collision frequencies ca n be expressed as multidimen- 
sional integrals over momenta of colliding particles. IShternin and YakovlevI (l2007l) 
calculated these integrals in the weak-screening approximation and described the re¬ 
sults by simple analytical formulae (referenced in Table |2]). In the case of strongly 
degenerate ultra-relativistic leptons, which is typical for neutron star cores, the latter 
authors obtained a much simpler expression, which can be written as 


Ke,p Ri20.8c(/?Fe,^/^)^- 


(17) 


The latter simplification, however, does not hold if the protons are superfluid. 


3 Superfluidity and superconductivity 


Soon after the development of the BCS theory (iBardeen et al. 


plains superconductivity by Cooper pairing of fermions (ICoonei 


1957). which ex- 


1956h . iBohr et al 


( 1958b argued that the same phenomenon of pairing is occurring inside nuclei (later 
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this suggestion was confirm ed experimentally). Mied^ j_959h extended the idea to 
the interior of neutron stars. Ginzburg and Kirzhnits formulated a number of 

important propositions concerning neutron superfluidity in the interior of neutron 
stars, the formation of Feynman-Onsager vortices, a critical superfluidity tempera¬ 
ture {Tc < 10^** K) and its dependence on the density (p ^ 10*^ -10*^ g cm^^), and 
discussed the influence of neutron superfluidity on he at ca pacity and theref ore on the 
thermal evolution of a neutron star. Bavm et HI ( 19691) and loinzburg ( 1970l) analyzed 
the consequences of neutron superfluidity and proton superconductivity: rotation of 
the superfluid component in the form of quantized vortices and splitting of the in¬ 
ternal stellar magnetic field into fluxoids (Sect. 14.3.41 ). Later many different authors 
considered various types of pairing of nucleons, hyperons, or quarks using different 
model potentials. 

Although we will not consider exotic models of neutron star cores, let us 
mention that superfluid i ty is possible in these models as well. For instance. 


Takatsuka and Tamagakil (1 19951) reviewed calculations of neutron and proton super¬ 


fluid gaps in pion condensed matter. Some authors have discussed superfluidity in 
quark matter (e.g^ Steiner et al. 20091). If h yperons are present, they can also be in a 
superfluid state ( Balberg and Barnealll998h. For a d e tailed recent review of superflu¬ 
idity in the interiors of neutron stars, see Page et aP ( 2014 ). 


3.1 Pairing types and critical temperatures 


The Cooper pairing appears as a result of the attraction of particles with the anti¬ 
parallel momenta,which is expected to occur, at low enough temperature, in any de¬ 
generate system of fermions in which there is an attrac tive interaction between par¬ 
ticles whose momenta p lie close to the Fermi surface (ICooDeiill956h . The strength 
of the interaction determines the critical temperature at which the pairing phase 
transition will occur. In a normal system the particle energy e varies smoothly when 
the momentum crosses the Fermi surface, while in the presence of pairing a discon¬ 
tinuity develops, with a forbidden energy zone having a minimum width of 2Apair at 
p = pp, which can be regarded as the binding energy of a Cooper pair. 

The BCS equations that describe symmetric nuclear matter in atomic nuclei and 
asymmetric neutron-rich matter in neutron stars have much in common but have also 
some differences. For instance, pairing in atomic nuclei takes place in the singlet 
state of a nucleon pair. In this case, the energy gap is isotropic, that is independent of 
the orientation of nucleon momenta. On the other hand, one can expect triplet-state 
pairing in the neutron-star matter, which leads to anisotropic gap. Singlet-state neu¬ 
tron superfluidity develops in the inner neutron star crust and disappears in the core, 
where an effective neutron-neutron singlet-state attraction transforms into repulsion. 
Triplet-state neutron superfluidity appears in the neutron-star core. Protons in the core 
can undergo the singlet-state pairing. 

The triplet pair states may have different projections nij of the total pair momen¬ 
tum onto the quantization axis: \mj\ = 0, 1, and 2. The actual (energetically favor¬ 
able) state may be a superposition of states with different my. Owing to uncertainties 
of microscopic theories this state is still unknown; it depends possibly on density and 
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Fig. 5 Critical temperatures of singlet neutron (left panel), triplet neutron (middle panel), and sin¬ 
glet proton (right panel) superfluidities in the inner crust and core of a neutron star, as func¬ 
tions of_gr^itational mass density p, for different superflu idity mod els, as marked near th e curves 
(see |Ho et |2015h: AO ^Amundsen and pstgaard [l985a|b|), AWP2 jAinsworth et alj Il989[), BCLL 
jBaldo et alJl992h. BEEHS iBaldo et alJl998h. BS iBaldo and Schulzdl2007h. CCDK jChen et all 199.31) . 

GIPSF lO aSoifi et alj|20d3). MSH~ lMareueron et al.ll20o3) ~SFB jSchwenk et al]|2003l) . TTav and TToa 

iTakatsukaanTldmagai^2004> . 


temperature. In simulations of neutron star cooling, one usually considers the triplet- 
state pairing with \mj\ =0 and 2, since the ir effects on the heat capacity and neutrino 
luminosity are qualitatively different (e.g.. I Yakovlev et alJll999bll200ll) . 


The critical temperature is very sensitive to the strength of the repulsive core 
of the nucleon-nucleon interaction. It is r elated to the superfluid ener gy gap by 


Tc = 0.5669zipair for the singlet gap (e.g., Lifshitz and Pitaevskh 20021 Sect. 40). 


For t he triplet gap, the situation is more complicated, because the gap is anisot ropic 


(e.g., Amundsen and 0stgaardllT985btlBaldo et al.l[l992l Yakovlev et al.l[iV99bl) . Ex¬ 
amples of the dependence of on gravitational mass density in the crust and core 
of a neutron star are sho wn in Fig. |5] Here, we employed the gap parametrization 
of Kaminker et al. ( 200lh with t he param e ter va lues and notations for different mod¬ 
els of superfluidity according to Ho et al.l ( 2015 ) togeth er with the p-depende nces of 
free-nucleon number densities n„ and rip from the fits dPotekhin et al.l 2013 ) for the 
BSk21 model of crust and core composition. Figure|5]demonstrates a large scatter of 
theoretical predictions, but also general features. We see that the ^So superfluidity 
of neutrons occurs mostly in the inner crust and the ^P 2 superfluidity mostly in the 
core. The critical temperatures of neutrons in the triplet states, 7’c„(^P2), and protons, 
7cp('So ), have usually a maximum at a supranuclear density p > po- Typical mag¬ 
nitudes of Tc vary from one model to another within a factor of a few. Neutron ^P 2 
superfluidity has, in general, much lower Tc than * So pairing of neutrons in the inner 
crust and protons in the core. 
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3.2 Superfluid effects on heat capacity 


Once a component x of the neutron star matter becomes superfluid, its specific heat 
Cv,x is strongly altered. When T reaches Te x, the critical temperature for the pairing 
phase transition, Cv,x jumps by a factor > 2. However, as T continues to decrease, 
the heat capacity becomes progressively suppressed. At T <C Ti the energy gap in the 
nucleon spectrum strongly reduces the heat capacity even compared to its value in the 
absence of pairing. These effects are implemented in numerical calculations through 
“control functions”/?c(7’/7c,x) as 


- JVC Ov.X 1 


(18) 


where denotes the value in the normal phase, Eq. (fTOl i. The control function de 


pends on the type of pairing. This dependence was studied bv iLevenfish and Yakovlev 
( 19941) . Analytical fitting formulae for in the npejX matter for the ma in types of 
superfluidity listed above are given by Eq. (18) of Yakovlev et al.l ( 1999bh i| 


Thr ee examples of the control functions, calculated according to I Yakovlev et al 


(Il999bl) (with the correction mentioned in footnote [T]l, are shown in the left panel 
of Eig. |6] One can notice that Cv,x nearly vanishes when T drops below ^ O.lTc x. 
Therefore, in the case of extensive pairing of baryons, the heat capacity of the core 
can be reduced to its leptonic part. This would result in a drastic reduction of the 
total specific heat, as already demonstrated by the heavy long-dashed line in Eig. [1] 
where we adopted MSH, TToa (assuming nij = 0), and BS superfluidity models for 
neutrons in the crust and core, and protons in the core, respectively, according to the 
notations in the caption to Eig. |5] 

Another example of the distribution of Cv among the various core constituents is 
shown in the right panel of Eig.|6l Here, we have adopted SEB, BEEHS (with mj = 0), 
and BCLL pairing gaps. The behavior of Cv as function of p proves to be qualitatively 
similar for different sets of superfluid gap models. In all cases this behavior strongly 
differs from that for unpaired nucleons, which is shown by thin lines for comparison. 


3.3 Superfluid effects on neutrino emission 


The enormous impact of pairing on the cooling comes directly from the appearance 
of the energy gap at the Eermi surface which leads to a suppression of all pro¬ 
cesses involving single particle excitations of the paired species. When T <C 7); the 
suppression is of the order of and hence dramatic. Its exact value depends 

on the details of the phase space involved in each specific process. In numerical cal¬ 
culations it is introduced as a control function. As well as for the heat capacity, for 
the neutrino emissivity one writes 


n _ E)(P“cesstype)^(0) 

^(pairing type) 5 


(19) 


^ In the latter paper, an accidental minus sign in front of the term (0.2846v)^ in the denominator of the 
fitting formula for Rc in the case of “type C” (^P 2 , \fnj \ — 2) superfluidity must be replaced by the plus 
sign (D.G. Yakovlev, personal communication). 
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stellar volume [10^® cm^] 
6 5 4 3 2 1 0 




Fig. 6 Left panel: Specific heat control functions for the ^So , (wty = 0), and ^P 2 i\>nj\ — 2) types 
of pairing listed in Sect. HI] The inset displays the same functions on a logarithmic scale. Right panel: 
Total and partial specific heats near the bottom of the crust and in the core of a neutron star at T = 10^ K 
as functions of density. The solid lines show the total Cy, and the other lines show the contributions of 
electrons (e~), neutrons (n), nuclei {N, in the crust), muons and protons {fi~ and p, in the core). Thin 
lines show results of a calculation with nucleons assumed to be unpaired, and thick lines take pairing into 
account. The top axis shows the volume contained inside a sphere with given p for a 1.4 Mq neutron star. 
The stellar stmcture and composition are adopted from the BSk21 model. 


where relates to the same process in the absence of pairing. These control func¬ 
tions (reduction factors) are available in the form of analytical fits, referenced in 
Tablell] 

The superfluidity not only reduces the emissivity of the usual neutrino reactions 
but also initiates a specific “pair breaking and formation” (PBF) neutrino emission 
mechanism. The superfluid or superconducting condensate is in thermal equilibrium 
with the single particle (“broken pairs”) excitations and there is continuous formation 
and breaking of Cooper pairs. The formation of a Cooper pair liberates energy which 
can b e taken away by a v — V pair ( Howers et al.lll976 ; Voskresensky and Senatorov 


1987h . This effect is most pronounced near the Fermi surface. When T falls below 


Tc, the neutrino emissivity produced by the Cooper pairing sharply increases. The 
PBF mechanism is sensitive to the model adopted for calculating the superfluid gaps 
in the baryon spectra: it is more important for lower gaps (weaker superfluid). Its 
emissivity is a sharp function of density and temperature. The main neutrino energy 
release takes place in the temperature interval between ~ Tc/S and T^. The control 
functions and the inten sity of the Cooper-pair neutrino emission are available as ana¬ 
lytical fits collected by Yakovlev et al.l (l200ll) (see references therein for the original 
der ivations), as indicated in Table [T] abo ve. 

Voskresensky and Senatorov d 1987b noticed that the PBF mechanism is sensi¬ 


tive to the in-medium renormalization of the nucleon weak-interaction vertex due to 
strong interactions (cf. Sect. 12,6.31 1. Later this ef fect has been reexami ned in many 
papers for different types of baryon pairing - see LeinsonI ( 20091 2010b for modern 
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results and a critical analysis of previous works. The net result is that the collective 
effects virtually kill down the PBF emission for the singlet pairing of nucleons, but 
leave this mechanism viable for the triplet pairing. Quantit atively, PBF emissivity es¬ 
timated without in-medium effects ( Yakovlev et aP 1999ah has to be multiplied by a 
small factor of [p^/m*cY \n the. 'Sq case, but by a moderate factor of ft! 0.19 in the 
^P 2 case. This result lies at the basis of the “minimal cooling scenario” and the expla¬ 
nation of the observed fast cooling of the neutron star in the Cassiopeia A supernova 
remnant (see Sect.|6]l. 

Superconductivity of protons may also induce another type of neutrino emission, 
electron-fluxoid scattering, in the presence of a strong magnetic field. It will be ad¬ 
dressed in Sect. 1431 


3.4 Superfluid effects on heat conduction 


The effects of nucleon superfluidity on the heat transp ort in neutron stars were dis¬ 
cussed qualitatively by Flowers and Itoh (197^, 1981- )■ The thermal co nductivity of 
electrons and muons was rec onsidered by iGnedin and YakoylevI (1199511 and later by 
Shternin and Yakoyleyl (1200711 . who obtained accurate analytical expressions yalid for 
a wide class of models of superfluid and non-superfluid matter. Baiko et al. ( 2001ah 
reanalyzed the thermal conduction by neutrons, utilizing some new deyelopments 
in the nucleon-nucleon interaction theory. The latter authors showed that the low- 
temperature behayior of the nucleon thermal conductiyity is yery sensitiye to the re¬ 
lation between critical temperatures of neutrons and protons. 

The lepton heat conduction in the core can also be affected by proton super- 
conductiyity, because superconductiyity modifies the transyerse polarization func- 
tion and screening functions in neutron-star matter. The se effects were stud ied by 
Shternin and Yakoyley ( 2007h . These authors, as well as Baiko et al. ( 2001 ah . man¬ 
aged to describe the effects of superfluidity by analytical functions, which facilitate 
their inclusion in simulations of neutron-star thermal eyolution (see Table|2|i. 

In the presence of neutron superfluidity, there may be another channel of heat 
transport, the so-called conyectiye counterflow of the normal component of matter 
with respect to the sup erfluid one. This mecha nism is known to be quite effectiye in 
superfluid helium (e.g., Tilley and Tilley! 1990li . but in the context of neutron stars the 
situation is unclear and has not been studied in sufficient detail. 

Heat can also be transported through the neutron star crust by collectiye in odes 
of superfluid neutron matter, called superfluid phonons ( Aguilera et al. 2009ll . At 
p ft 10*^ — lO'^ g cm^^ the conductiyity due to superfluid phonons was estimated 
to be significantly larger than that due to lattice phonons and comparable to electron 
conductiyity when T ft 10^ K. The authors found that this mode of heat conduction 
could limit the anisotropy of temperature distribution at the surface of highly magne¬ 
tized neutron stars. Howeyer, new stu dies of the low - energy collectiye ex citations 
in the inner crust of the neutron star ( Charnel 2012 : Charnel et al. 2013h . includ¬ 
ing neutron band structure effects, show that there is a strong mixing between the 
Bogoliuboy-Anderson bosons of the neutron superfluid and the longitudinal crystal 
lattice phonons. In addition, the speed of the transyerse shear mode is greatly re- 
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duced as a large fraction of superfluid neutrons are entrained by nuclei. This results 
in an increased specific heat of the inner crust, but also in a decrease of the electron 
thermal conductivity. On the other hand, the entrainment of the unbound neutrons 
decreases the density of conduction neutrons, i.e., neutrons that are effectively free. 
The density of the conduction neu trons can be much smaller than the total density of 
unbound neutrons ( Chamelll2012 ). which results in a decrease of the neutron thermal 
conductivity. 


4 The effects of strong magnetic fields 

4.1 Magnetic-field parameters 


Convenient dimensionless parameters that characterize the magnetic field in a plasma 
are the ratios of the electron cyclotron energy hcOc to the atomic unit of energy, elec¬ 
tron rest energy, and temperature: 


7 = 


h^B 


B 


b = 


heB 


^^ = ^ = 134 . 34 ^. 
Bqed ’ r 76 


( 20 ) 


Here, (Oc = eBjineC is the electron cyclotron frequency, Bq = 2.3505 x 10‘^ G is the 
atomic unit of ma gnetic field, ZloFn = 4.414x lO'^Gis the critical field in Quantum 
Electrodynamics ischwing^ 1988 ). and B \2 = B/IO'^ G. 

Motion of charged particles in a magnetic field is quantized in discrete Lan¬ 
dau levels. In the non-relativistic theory, the energy of an electron in a magnetic 
field equals Nh(Oc -\-mep\l2, where is the momentum component along B, N = 
wl + 5 T 5 characterizes a Landau level, the term T 5 is the spin projection on 
the field, and is the non-negative integer Landau number related to the quanti- 
zation of the kinetic moti on transverse to the field. In the relativistic theory (e.g., 
Sokolov and Terno^ll986h . the kinetic energy £ of an electron at the Landau level N 
depends on its longitudinal momentum as 


^n[Pz) = c {m^c^ + IhcOcirieN + pfj — mgC^. 


( 21 ) 


The levels ^ 1 are double-degenerate with respect to the spin projection s. Their 
splitting de due to the anomalous mag netic moment of the electron is negligible , 
because it is much smaller than h(Oc (e.g., Schwing^ll988 : Suh and Mathewsll2001 ): 


~ ^ Of ( hcOc at <C 1 , 

^ \ nteC^ [ln/7 — 1.584]^ 1, 


( 22 ) 


where ctf is the fine structure constant. 

The Landau quantization becomes important when the electron cyclotron energy 
hcOc is at least comparable to both the electron Lermi energy £f and temperature T. If 
hcOc is appreciably larger than both £f and T, then the electrons reside on the ground 
Landau level, and the field is called strongly quantizing. The condition hcOc > T is 
equivalent to > 1. The condition TicOc > Ep e translates into p < ps, where 

Pb ^ 7045 Y-^ B%^ gcm-\ 


( 23 ) 
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In the opposite limit, where either <C 1 or p ps, the field can be considered as 
nonquantizing. 

For the ions, the cyclotron energy is h(Oa = Z {me/mi)h(Oc, and the Landau quan¬ 
tization is important when the parameter 


Ci = hco,i/T = 0.0737 {ZlA)BnlT(. 


(24) 


is not small. The energy spectrum of an ion essentially differs from Eq. (EB because 
of the non-negligible anomalous magnetic moments. In the non-relativistic theory, 
the energy of an ion equals e = (wl + + tn\p^/'^ + where hl is the 

ion Landau number, is the longitudinal momentum, gi is the g-factor (gi = 2 in 
the Dirac theory, but, e.g., gi = 5.5857 for the protons), and si is the integer quantum 
number corresponding to the spin projection on B in units of hjl. If the ions are 
relativistic, the situation is much more com plicated. For baryons w ith spin ^ (e.g., 
protons) the energy spectrum was derived by Broderick et al. (l2000[) . 


4.2 Magnetic field effects on the equation of state and heat capacity 
4.2.1 Magnetized core 


A magnetic field can affect the thermodynamics of the Coulomb plasmas, if the Lan¬ 
dau quantization is important, i.e., under the conditions that are quantified in Sect. 14. II 
In particular, Eq. (l23t can be recast into 

B > (3.8 X 10^^ G) {Yen^/im-^f/^. (25) 


We have n^^ 0.1 fm^^ near the crust-core interface, and Yg is typically several per¬ 
cent throughout the core. Therefore, the electron component of pressure in the core 
might be affected by the fields B > 10*^ G. 

One can easily generalize Eq. (|25]) for other fermions (p-mesons, nucleons) in the 
ideal-gas model. In this case. Ye should be replaced by the number of given particles 
per baryon, and the right-hand side should be multiplied by = 206.77 for 

muons and ^10^ (of the order of nucleon-to-electron mass and electron-to-nucleon 
magnetic moment ratios) for protons and neutrons. Accordingly, the partial pressures 
of muons and nucleons in the core cannot be affected by any realistic {B <a few x 
10^ ^ G) magnetic field. 


Broderick et al. ( 20001) developed elaborated models of matter in ultra-magnetized 


cores of neutron stars. They considered not only the ideal npep gas, but also interact¬ 
ing matter in the framework of the relativistic mean field (RMF) model. The magnetic 
field affects their EoS at B > 10*^ G. As follows both from the estimates based on 
the virial the orem (Lai and Shapird 199j) and from numerical hydrodynamic simu- 
' M 


lations (e.g.. iFrieben and Rezzollall2012L and references therein), this field is close 
to the upper limit on B for dynamically stable stellar configurations. The effect is 
even smaller when the magnetization of matter is included consistently in the EoS 
(IChatteriee et al.ll20L5h . Therefore, it is unlikely that a magnetic modification of the 
EoS could be important in the cores of neutron stars. 
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Fig. 7 Left panel: Normalized pressure Pln\T\ entropy S and heat capacity Cy per one ion, and logarithmic 
derivatives of pressure over density and temperature, Xp and Xt^ foi" ^ fully-ionized non-magnetic (dashed 
lines) and magnetized {B — 10^^ G, solid lines) iron plasma dXT — 10^ K. The vertical dotted lines mark 
the densities at which the electron Fermi temperature e quals T without (7) or with (2 ) the magnetic field, 
p — pg (i), rc—PjYi (4), and Tp — T (5). (Figure 6 from fpotekhin and Chabriejl20131 reproduced with the 
permission of ©ESO.) Right panel: Normalized thermal phonon contribution to the reduced heat capacity 
as a function of logio(7’ /Tp) at different values of the ratio hcOd/Tp, marked near the curves. 


4.2.2 Magnetized crust and ocean 


AtB> 10^® G, nuclear shell energies become comparable with the proton cyclotron 
energy. Thus the interaction of nucleon magnetic moments and proton orbital mo¬ 
ments with magnetic held may cause appreciable modihcations of nuclear shell en- 
ergies. These mo d ihcati ons and their consequences for magnetars were studied by 
Kondratyev et al.l (1200 ih . who found large changes in the nuclear magic numbers 
under the inhuence of such magnetic helds. This effect may alter signihcantly the 
eq uilibrium chemical com position of a magnetar crust. 

Muzikar et al. ( 1980h calculated the triplet-state neutron pairing in magnetized 
neutron-star cores. According to these calculations, magnetic helds B > lO'® G make 
the superhuidity with nodes at the Fermi surface energetically preferable to the usual 
superfluidity without nodes. Accordingly, the superhuid reduction factors for the heat 
capacity and neutrino emissivity (the control functions) may be different in ultra- 
str ong helds. 


Charnel et alJ (120121) studied the impact of superstrong magnetic helds on the 


composition and EoS of the neutron star crust. In particular, they found that the 
neutron-drip pressure increases almost linearly by 40% from its zero-held value in 
the interval lO'® G < B < 5 x lO'® G. With further increase of the held strength, the 
drip pressure becomes directly proportional to B. 

Thus the ultra-strong helds B > 10^® G can affect various aspects of the physics 
of the inner crust in quite non-trivial way. Hereafter we will consider only helds 
B < 10*® G. They can be quantizing in the outer crust of a neutron star, but not in the 
inner crust or the core. Analytical htting formulae for the thermodynamic functions 
of the electron-ion plasmas in such helds, as well as a computer code that imple- 
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ments these fitfl were published by Potekhin and Chabrieij ( 2013 ). Such fields affect 
the electron part of thermodynamic functions in the outer envelopes only, as illus¬ 
trated in the left panel of Fig. |7]in the case of fully-ionized iron at F = 10^ K and 
B — 10*^ G (for illustration, the density range is extended to p < 10^ neglecting 
the bound states that can be important in this p-T domain). We plot the principal 
thermodynamic quantities normalized per one ion as functions of density. For com¬ 
parison we also show them in the absence of quantizing magnetic field. The vertical 
dotted lines marked by numbers separate different characteristic domains, consecu¬ 
tively entered with increasing density: onset of electron degeneracy at B = 0 and at 
B = 10*^ G, population of excited Landau levels (p = p^), melting point with for¬ 
mation of a classical Coulomb crystal (I'm = T), and onset of the quantum effects in 
the crystal (IJ, = T). The gradually decreasing oscillations correspond to consecutive 
filling of the electron Landau levels. The magnetic field B = 10*^ G does not affect 
the ion contributions at this T. 

The contributions of ions to the thermodynamic functions are affected by the 
magnetic field if the parameter i^i, defined by Eq. (l24l l. is large. This may occur in 
a superstrong field of a magnetar. The right panel of Fig. [T] illustrates the effects 
of a s uperstrong field on Cv,i for the model of a harmonic Coulomb crystal ( Baikol 


20091) . Flere we plot the thermal phonon contribution to the heat capacity of the bcc 


Coulomb lattice calculated as the derivative Cy = TdS/dT of the fit to the phonon 
entropy S given by Eq. (77) of [Potekhin and Chabrien (120131) . This approximation is 
more accurate for the heat capacity than the alternative approximation that provides 
exact fulfillment of the Maxwell relations (Eq. 80 of the same paper). The three steps 
on the curves in the right panel of Eig.|7]correspond to contributions of three branches 
of the phonon spectrum, which are affected differently by the quantizing magnetic 
field. 


4.3 Magnetic field effects on neutrino emission 
4.3.1 Magnetic Durca process 


We have mentioned in Sect. 12.6.21 that the Durca reaction is the most efficient neu¬ 
trino emission process, but it can only operate above a certain t hreshold density in 
the cen tral parts of the cores of sufficiently massive neutron stars. iLeinson and Perez 
(Il998n noted that a superstrong magnetic field can substanti ally weaken this re¬ 
quirem ent. An accurate study of this effect was performed by iBaiko and Yakovlev 
( 19991) . They showed that the border between the open and closed Durca regimes 
is smeared out over some B-dependent scale and described this smearing by sim¬ 
ple formulae. In practice this effect should be very important for neutron stars with 
B > 10^® G. At less extreme fields (lO^'* G < B < 10*® G) it is important for neutron 
stars whose mass happens to be c lose (within a few percent) to the Durca threshold 


mass. 


Baiko and Yakovlev ( 1999h also showed that a strong magnetic field has a non¬ 


trivial effect (oscillations of the reaction rate) in the permitted domain of the Durca 
reaction, but the latter effect, albeit interesting, appears to be unimportant. 


Also available at http: //www. iof fe.ru/astro/EIP/. 
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Fig. 8 Neutrino emissivity in a magnetized crust from the synchrotron processes for two, uniform, mag¬ 
netic field strengths of 10^^ G (left panel) and 10^^ G (right panel). The contour lines are labeled by the 
value of logio[2v/(ergcm“^ s“^)]. Regions where this process dominates over the ones shown in Fig. [3] 
are lightly shadowed (in yellow) and regions where it dominates by more than a factor of 10 are darkly 
shadowed (in orange). The two dotted lines show the dominance transitions between the three processes 
presented in Fig. [3 (Also indicated is the ion melting curve, dashed line.) 


4.3.2 Pair annihilation 


The e pair annih ilati on process in strong magneti c fields 


s studied by 

Kaminker et al.l (Il992h and iKaminker and Yakovlevl(ll994l) . In a hot, non-degenerate 
plasma {T > lO*** K) only ultra-strong magnetic fields B > 10*® G can significantly 
affect the neutrino emissivity. Such fields can be quantizing in the p — T domain 
where the pair emission dominates (see Fig.O. They amplify Qpair by increasing the 
number densities of electrons and positrons via very strong quantization of their mo¬ 
tion. Lower fields may also influence Qpair but less significantly. A field B ^ 10*^ G 
may quantize the motion of positrons at T < 10® K and increase the positron number 
density. In this way the presence of a strong magnetic fields greatly enhances Qpair in 
a not too hot plasma. However, this enhancement usually takes place where the pair 
annihilation emissivity is much lower than the contribution from other neutrino reac¬ 
tions, and therefore it is unimportant for studies of neutron-star thermal evolution. 


4.3.3 Synchrotron radiation 


A relativistic electron propagating in the magnetic field can emit neutrinos because 
of its rotation around the magnetic field lines. This process is quite analogous to the 
usual synchrotron emission of photons. The calculation of the corresponding neutrino 
e missivity, Qsvn, is s imilar t o that of the pa i r annih ilatio n process. It was studied, e.g., 
by iKaminker et al. ( 1992h . Vidaurre et al. ( 19951) . and Bezchastnov et al.l (1997). In 
Fig. Owe show the plot of Qsyn on the p — T plane for two field strengths typical for 
magnetars, B = lO*'* G and 10*® G. It is clear from this plot that the the synchrotron 
process can be dominant in the crust of magnetars in a large temperature range. 
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4.3.4 Electron-fluxoid scattering 

The internal stellar magnetic field can be confined in the crust or be distributed over 
the entire star. In the latter case, a transition to a superconducting state in the course 
of stellar cooling is accompanied by a dramatic change in the spatial structure of 
the magnetic field. Initially homogeneous field splits into an ensemble of Abrikosov 
fluxoids - quantized magnetic flux tubes, which contain a superstrong magnetic field, 
embedded in the field-free superconducting medium. Neutrino synchrotron radiation 
is then modified and may be treated as neutrino parr emiss ion due to scatt e ring o f 
electrons on the fluxoids. This mechanism was studied by iKaminker et al.l (119971) . 
who obtained an analytical fit to the corresponding neutrino emissivity (referenced 
in Table [T]). The concentration of the field within the fluxoids amplifies the neutrino 
emissivity, compared to the usual synchrotron regime, when temperatures drops be¬ 
low the critical one for the protons, T^p. As long as T is not much lower than T^p, 
the Cooper pairing mechanism remains much more powerful (unless B ^ 10*® G, 
in which case the electron-fluxoid scattering may be more powerful at any T). At 
T Tcp, the electron-fluxoid scattering becomes the dominant neutrino emission 
mechanism for the neutron stars with strong and superstrong fields {B > 10*^ G). 


4.4 Magnetic field effects on heat conduction 
4.4.1 Photon heat conduction 


The thermal conductivity K is related to the opacity k by the equation 


K = 


3p>c ’ 


(26) 


where cjsb is the Stefan-Boltzmann constant. The spectral radiative opacities for two 
normal polarization modes in strongly magnetized neutron-star photospheres are re¬ 
viewed in Paper I. These opacities, >C(pj{9b), where j = 1,2 marks the extraordinary 
and ordinary polarization modes, depend on the angle 9b between the wave vector 
and magnetic field. In the diffusion approximation, they combine into the effective 
opacities for the transport along p and across (x^ j) magnetic field according to 


( 4 .;) ' 




Jo 


2 cos^ 9 b 1 sinP^d^B 
singes j Xa,j{9B) ' 


(27) 


The effective opacity for energy transport at angle 0 to B in each polarization mode 
is given by 1/= cos^ 9 + sin^ 9 /where x is the Rosseland mean of >C(p, 


Jo 


u{z) ^ ^ 15 

-dz, u{z) = - 

xta.j 




hco 


z = 


4a:4 (e^- 1)2’ 


(28) 
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For fully ionized plasmas, the radiati ve opacities are contributed fr om the free- 
free absorption and Thomson scattering. Silant’ev and Yakovlev! ( ToSOh studied the 
Rosseland opacities for a non-polarized radiation in magnetized fully ionized plasmas 


= 


-I- 


(29) 


using the Born approximation for the free-free contribution. iPotekhin and Yakovlev 
( 2001 ) obtained simple analytical fits for including a correction to the Born 


approximation, as functions of p,T,Z, A, and the magnetic-field parameter defined 

by Eq. (l20l) . Asymptotically, at —>■ 

At finite but large the radiative opacities of fully ionized matter are strongly 
reduced. The reduction is ~ 10 times stronger for the Thomson scattering than for the 
free-free absorption. In deep, strongly magnetized photospheric layers the Thomson 
scattering do mnates only if Tf, > 10p^(^ > 10 Bj 2 ^; otherwise the free-free absorp¬ 
tion prevails (IPotekhin and Yakovlevl200lh . 

The outermost envelopes of neutron stars can be incompletely ionized in the cases 
of large Z or B. The presence of bound species can strongly affect the radiative opac¬ 
ities and the spectrum of emitted radiation, as discussed in Paper I. However, the 
layers that are responsible for the heat flux from the interior of the neutron star to the 
surface, as a rule, lie at sufficiently large densities, where the plasma is fully ionized 
by pressure. Therefore the bound species are usually unimportant for the effective 
surface temperature of a neutron star. 


4.4.2 Electron heat conduction 


A non-quantizing magnetic field does not affect thermodynamic functions of the 
plasma. However, it does affect the electron heat conduction, if the Hall magneti¬ 
zation parameter 


CUgT 


1760 


Bn 

■\/l -\-x^ 


T 

10-1® s 


(30) 


is not small. Here, cOg = Wc/v/l +xj is the electron gyrofrequency, and T is the 
effective relaxation time. In a degenerate Coulomb plasma with a non-quantizing 
magnetic field, the main contribution is given by the electron-ion scatterin g according 
to Eq. (IA.3b . This regime has been studied by Yakovlev and Uminl ( 1980ll . 

Electric and thermal currents induced in a magnetized plasma under the effect 
of an electric field E, a weak gradient V/r of the electron chemical potential, and a 
weak temperature gra dient VT can be decomposed into conduction and magnetiza¬ 
tion components (e.g., Hernauis^ 1984 ). The latter ones relate to surface effects and 
must be subtracted. Let and jj be the conduction components of the electric and 
thermal current densities. They can be written as 


j, = d-E*-a-VT, jT=P-E*-KVT, (31) 


where E* = E + Vp./e is the electrochemical field. The symbols a, a, j3, and k 
denote second-rank tensors (<7 is the conductivity tensor) which reduce to scalars at 
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B — 0. Equations dSTT i can be rewritten ; 


E* =Rj,-Q-VT, j^ = -TQj,-kVT, 


(32) 


where R = Q = —R ■ a, and K'=K'+7’d:-2are the tensors of specific resistance, 
thermopower, and thermal conductivity, respectively. 

Electron heat and charge transport controlled b y electron-ion collisions in quan ¬ 
tizing magnetic f i elds of neutron sta r s was stud i ed by Kaminker and YakovlevI (Il98lh ; 
Yakovlev ( 1984 ): Hernauis^ t 1984 ): Potekhin ( 1996' fl999h . The components of ten¬ 
sors ff, a, and k can be expressed as (|Potekhinll99^ 


an > = 


I 


e{ll-s)/T 


^sje) 




1 


( 9 e e (® -|-1 


de, (33) 


where 


^(e) = 


Ar=0 


Ymax 

(2- 5Ar,o)|p2|, 


(34) 


Mnax is the maximum Landau number for a given electron energy e, and \p^\ depends 
on e and N according to Eq. (12111 . In a non-quantizing magnetic field, i.e., at A^max ^ 
1, the sum can be replaced by the integral, which gives ./fe(e) = (p/^)^/37r^, where 
p is the momentum that corresponds to the energy e. The functions Tij{e) play role of 
relaxation times for the components of tensors ff, a, and k, determined by electron 
scattering. In general, they differ from the mean free time Tei(e) = 1 /Vei{s) between 
scattering events for an electron with energy e. Because of the symmetry properties 
of the tensors ff, a, and k, in the coordinate frame with z axis directed along B, 
there are only three different non-zero components of Tif related to longitudinal 
currents, = Xyy related to transverse currents, and Xxy = —Tyx related to the Hall 
currents. 

In a quantizing magnetic field, there are two different effective mean-free times 
T|[ (e) and Tx (e), corresponding to electron transport parallel and perp endicular to B. 
In this case, the classical expressions (e.g., Yakovlev and IJminll 19801) are recovered: 


Tzz = Til 


To: — 


Tx 


l + (®gTx)2 


Tyx — 


CUgTi 


l + (®gTx)2 


(35) 


It is convenient to keep usin g Eg. (IA.3I ) for T[[ and Tx, but with different Coulomb 
logarithms A|| (e) and A i (el. lPotekhinl (Il999h calculated these Coulomb logarithms 
and fitted them by analytic expressions. Their Fortran implementation is available at 
http: //www. iof f e. ru/astro/conduct/. In the limit of non-quantizing magnetic 
field, T|| = Tx = Tgi(£) is given by Eq. (IA.3I) with e = jX. 

When the electrons are strongly degenerate, the derivative in the square brackets 
in Eq. (l33]) is sharply peaked. Then Eq. (l33t gives 


2 2 

e c He 




n^T 


( 36 ) 
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Fig. 9 Electron thermal conductivities along (upper curves) and across (lower curves) magnetic field B = 
10*^ G (left panel) and lO'^ G (right panel) as functions of mass density at temperatures T = 10^ K (solid 
lines) and 10* K (dot-dashed lines). For comparison, the non-magnetic thermal conductivities are shown 
by dotted lines. 


The latter relation is the Wiedemann-Franz law generalized to the magnetic case. On 
the other hand, Eq. (l3Tt satisfactorily describes the conductivities in general, includ¬ 
ing the opposite case of weakly degenerate electrons. 

Figure|3illustrates the p -dependence of the thermal conductivities along (K'||) and 
across {k±) the magnetic field. The first, most significant peak at JCn is related to the 
filling of the first Landau level by the electrons atp^ps- The other peaks correspond 
to consecutive filling of higher Landau levels. 


5 Thermal structure of neutron stars 

5.1 Blanketing envelopes 

The very different thermal relaxation timescales of the envelope and the crust of a 
neutron stars makes computationally expensive to perform cooling simulations in a 
numerical grid that comprises both regions. Radiative equilibrium is established in the 
low-density region much faster than the crust evolves, so that the envelope reaches 
a stationary state on shorter timescales. Thus, the usual approach is to use results of 
stationary envelope models to obtain a relation between the photon flux Fph radiated 
from the surface and the flux Ft, and temperature T\, at the crust/envelope boundary, 
p = Pb- This relation supplements the evolution equations for the interior [Eq. 0] as 
an outer boundary condition. 

The boundary density Pb is chosen as a trade-off between two requirements; first, 
that the thermal relaxation time of the layer with p < Pb is short compared to the 
characteristic variability time of the studied thermal radiation, which favors smaller 
Pb, and second, that T does not strongly vary at p > pb, which favors larger pb. For 
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weakly magnetized, isolate d cooling neutron stars, pb is usually set at 10“' g cm 


10 —3 

^_ _ _ _ _ ivL tiL iu gem 

( Gudmundsson et ^Il983 ). but in general it varies from 10^ g c m~^ for neutron stars 
with rapid variations of thermal emission ( Shternin et al.ll2007h to Pdup for relatively 
hot and strongly magnetized neutron stars ( Potekhin et alj|2003 ). 

At every Tb, Fph or, equivalently, the effective surface temperature T^, depends 
on the properties of the heat-blanketing envelopes. In the absence of neutrino energy 
losses in the envelope (that is the case for most cooling neutron stars, except for the 
hottest ones), the flux Fph at the surface is equal to the flux Ft, at the inner boundary of 
the blanketing envelope. Then it is sufficient to know the 7), - relation for cooling 
sim ulations. _ _ 

Gudmundsson et ^ ( 1983 ) carried out a comprehensive study of the thermal 


structure of the non-magnetized blanketing envelopes composed of iron, using 
the best physics input available at that time. They considered the envelopes with 
logjQ7’s[K]> 5.25 (there were no reliable calculation of the thermal conductivities 
for lower temperatures) and fitted the numerical solutions by a remarkably simple 
formula 

Tb = 1.288 X 10^ K, (37) 

where Ff, = Ts/lO^ K . An a nalytical derivation of a similar expression was given by 
Ventura and PotekhinI (1200 ih . A more accurate but less simple fit was constructed by 
Potekhin et al. ( 1997 ). 

The Th - Ts relation is mainly regulate d by the thermal conductivity in the “sen¬ 


sitivity strip” (iGudmundsson et al.lll983h that plays the role of a “bottleneck” for 


the heat leakage. Its position lies around the line where jq- = jq, (as a rule, around 
p - 10^-10^ g cm ^ for B — 0) and depends on the stellar structure, the boundary 
temperature 7b, the magnetic held B in the vicinity of the given surface point, and the 
chemical composition of the envelope. Since the magnetic held hampers heat trans¬ 
port across B, the depth of the sensitivity strip can be different at different places of a 
star with a strong magnetic hel d: it lies deeper at the places where the magnetic held 
lines are parallel to the surface ( Ventura and Potekhinll2001 ). 

The blanketing envelopes are more transparent to the heat flux, i f they are com 


posed of light chemical elements. This eff ect was studied in det ail by (Potekhin et al 


( 1997 ) for non-magnetic envelopes and bv Potekhin et al.l ( 2003 ) for strongly magne¬ 
tized envelopes. The effect is related to the Z-dependence of the collision frequencies 
Ve\- The higher is Z, the larger is Vei and the lower is the conductivity. A temperature 
variation by a factor of 30 can change the thermal conductivity of iron plasma less 
than altering the chemical composition from Fe to He at a fixed T. This effect has 
important consequences for the relationship between the surface and internal tem¬ 
peratures of neutron stars. For example, combined effects of strong magnetic fields 
and light-element composition simplify the interpretation of magnetars; these effects 
allow one to interpr et observations assuming less extreme (therefore, more rea listic) 
heating in the crust ( Kaminker et al.ll2009t Pons et aDl2009l: Vigano et al.ll2013 ). 

The envelope is thin (its depth Zb 100 m, if Pb = 10**^ g cm^^) and contains a 
tiny fraction of the neutron-star mass (~ 10^^, if Pb = 10*** g cm^^). Therefore one 
can neglect the variation of the gravitational acceleration in this layer. Neglecting 
also the non-uniformity of the energy flux through the envelope due to the neutrino 
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emission (which is small, if the neutron star is not too hot, as we discuss below) and 
the variation of the temperature over the surface (which varies on larger length 
scales tha n Zh), one can obtain, instead of Eg. (B, the much simpler thermal structure 
equation (iGudmundsson et al.lll983l : IVan Ripeiill988h 


dlnr _ 3 Pk 
dlnP “ 16 


(38) 


wherex: is the total opacity, related to the conductivity K via Eq. (l26l l. 

The assumption of a constant flux, however, breaks down for magnetars, most of 
which have atypically high effective temperatures. In this case one should solve the 
complete set of equations, taking neutrino emission and heat sources into account. 
The neutrino emissi on from the crust lim its the effective surface temperature of a 
cooling neutron star (iPotekhin et al.ll2007l : cf. Fig. [TT]below). This ve ry same effect 
is wh at limits the maximum flux in the few days of a magnetar outburst (IPons and Rea 
20121). In addition, for magnetars one must go beyond the plane-parallel approxima¬ 
tion (see Sect. 15.3b 


5.2 The effects of strong magnetic fields 


As seen from Eqs. (l30b and (IA.3b . the Hall magnetization parameter is large in the 
outer neutron-star envelope at B > 10^* G. Moreover, the magnetic held can be 
strongly quantizing in the outermost part of the envelope. In this case the magnetic 
held can greatly affect the heat conduction and the thermal structure. 

Figure [To] shows example s of the temperature profiles in the envelopes. The left 
panel is a recast of Fig. 8 from Potekhin and Chabriei ( 2013h . Here we show a profile 
of an ultra-magnetized neutron star, with B = 10*^ G, and with relatively high surface 
temperature, log[o7’s(^r) = 6.5, which is similar to the values evaluated for some 
magnetars. In this case, thermal photons are radiated from a solid surface, with high 
mass density p —2x 10^ g cm^^ just below the surface. The temperature quickly 
grows at the solid surface and reaches the melting point at the depth z Ri 7 cm. Thus, 
at the given conditions, the liquid ocean of a magnetar turns out to be covered by 
a thin layer of “ice” (solid substance). We treat the solid crust as immobile, but the 
liquid layer below the “ice” is convective up to the depth z ^ 1 m. The change of 
the heat-transport mechanism from conduction to convection causes the break of the 
temperature profile at the melting point in Fig. [T0| We underline that this treatment 
is only an approximation. In reality, the superadiabatic growth of temperature can 
lead to a hydrostatic instability of the shell of “ice” and eventually t o its c racking and 
fragmentation into turning-up “ice floes”. Potekhin and Chabrieil ( 20131) speculated 
that such eyents may result in yariations of thermal luminosity of magnetars. The 
temperature profile flattens with density increase, and the Coulomb plasma freezes 
again at the interface between the layers of ®®Ni and ^®Kr at p = 1.5 x 10® g cm^^. 

For comparison, we also show the thermal profile without the magnetic field. It 
is smooth. There is neither magnetic condensation nor conyection. In this case, the 
spectrum is formed in the gaseous atmosphere at much lower density beyond the 
frame of the left panel. 
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Fig. 10 Thermal structure of blanketing envelopes with different magnetic helds. Left panel: temperature 
profile (solid line in the left top panel) for an envelope with ground-state composition, with Z values 
shown in the bottom left panel, for a neutron star with surface gravity g — 1.6 x 10^"^ cm s“^, magnetic 
field B — 10^^ G directed at 45° to the surface, and internal temperature Ti, — 6.7 X 10^ K, which yields 
the effective surface temperature Ts = 3.16 x 10^ K. The dot-dashed line is the melting line. The asterisks 
confine the part of the profile where heat is carried by convection. For comparison, the non-magnetic 
profile with the same is shown by dotted line. Right panel: temperature profiles for carbon blanketing 
envelopes for a neutron star with g — 1.4 x 10^^ cm s“^, = 10^ K, and magnetic fields 5 = 0 (dotted 

line), 10^^ G (short dashes), 10^^ G (long dashes), 10^"^ G (alternating short and long dashes), and 10^^ G 
(solid line). The vertical segment of the dot-dashed melting line corresponds to the maxim um density for 
carbon, according to the thermonuclear stability criterion of iPotekhin and Chabrien |20l3), where cai'bon 
gives way to heavier chemical elements, which form a crystal. The heavy dots mark the position of the 
radiative surface, where T equals the effective surface temperature T^. 


In the right panel of Fig. [TO] we compare temperature prohles for a neutron star 
with internal temperature 10^ K and heat blanketing envelopes made of carbon, en¬ 
dowed with different magnetic helds. For the held strengths up to 10*^ G, the radi¬ 
ation is formed in the gaseous atmosphere, whose density gradually becomes larger 
with increasing magnetic helds, due to the reduction of the effective opacities dis¬ 
cussed in Sect. 14.4.11 The temperature prohles are rather smooth. The blanketing 
envelopes are liquid at this temperature. At the largest held strength B = 10^^ G, 
however, the situation is qualitatively different. As well as in the case of the hotter 
ultra-magnetized ground-state envelope in the left panel, the heat is radiated from 
the condensed solid surface. Below the surface, at density 10^ g cm^^, the tempera¬ 
ture quickly grows, which causes melting of the Coulomb crystal with formation of 
a Coulomb liquid beneath the solid surface. With further density increase, the prohle 
suffers a break at ps « 4.5 x 10^ g cm^^ [Eq. (l2Tt l. where the electrons start to pop¬ 
ulate the hrst excited Landau level, which is associated with the peak of the thermal 
conductivity around ps (cf. Fig. 13. 

As we have seen in Sect. 14.41 the conduction is strongly anisotropic in these con¬ 
ditions. Therefore the effective local surface temperature is non-uniform and de¬ 
pends on the magnetic held geometry. Figure [TT] shows examples of the relations 
between and 7b deep in the crust for the magnetic helds B = 10^^ G and 10^^ G 
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perpendicula r and parallel to the ra dial direction. The relations obtained in the ID ap¬ 
proximation ( Potekhin et al.ll2()()7 ) with and without allowance for neutrino emission 
are plotted by the solid and dotted lines, respectively. We see that at Ti, < 10* K the 
neutrino emission does not affect Tg. At higher T\,> 10^ K, in contrast, this emission 
is crucial; if Qy — 0, then continues to grow up with increasing Tj,, whereas with 
realistic Qy the surface temperature tends to a constant limit, which depends on B. In 
most cases this limit is reached when Ti, ^ 10® K. 

Since the distribution of Ts over the neutron-star surface is non-uniform in strong 
magnetic fields, it is convenient to introduce the overall effective temperature of the 
star, Teff, defined by 


4;r(JsB^^7rff = Uh = 


^ph • 


j TphdZ = ctsb j 


tUe, 


(39) 


where is the local flux density and dE is the surface element. The quantities T^, 
Teff, and Lpi, refer to a local reference frame at the neutron-sta r surface. The redshifted 
(“apparent”) quantities as detected by a distant observer are llh orn 3l977h : 


R-^R/^l-rjR, TSi = T,n^l-rjR, l;^ = {I - rjR)]^^. (40) 

The effects of quantizing m agnetic fields on t he thermal structure of n eutron-star 
envelo pes were firs t studi ed by IHernauislI (Il985h and so mewhat later by Van Ripei 
( 1988 ) and Schaaj ( 1990ll . using the ID approximation. Van Ripei ( 19881) consid- 
ered a neutron star with a constant radial magnetic field. In this model, the quantum 
enhancement of conductivity at p near ps, seen in Fig. |3 results in an overall en- 
hancemen t of the neutron-star photon luminosity Lph at a fixed T\,. Consequently, 


Van Ripe^ (1991) found a stro ng effect of the magnet i c field B ^ 10** G on the 


neutron-star cooling. However, Shibanov and Yakovlev ( 19961) showed that, for the 
dipole field distribution, the effects of suppression of the heat conduction across B 
at the loci of nearly tangential field can compensate or even overpower the effect of 
the conductivity increase near the norm al direction of th e fie ld lines . Th^ conclu¬ 
sion confirmed the earlier conjectures of Hernauisl ( 1985 ) and Schaal ( 1990l) . In the 
2000s, detailed stu dies of the Tj, - F relation in str ong magnetic fields were performed 
for iron envelopes ( Potekhin and Yakovlevll200ll) and accreted envelopes composed 
of light elements ( Potekhin et al.l 2003 ). as wel l as for the large-sca le (dipole) and 
small-scale (stochastic) surfa ce magnetic fields (Potekhin et al.ll2005h . These studies 
confirmed the conclusions of Shibanov and Yakovlev (Il9% ). but showed that in su¬ 
perstrong fields B > 10*^ G the quantum enhancement of the conductivity and the 
corresponding increase of at the places where B is nearly radial overpowers the 
decrease in the regions of nearly tangential field lines, so that Tj-ff at a given in¬ 
creases. However, this may not be the case in the configurations where the field is 
nearly tangential over a signif icant portion of the stellar surface as, e.g., in the case 
of a superstrong toroidal field ( Perez-Azorin et allbOOb : Page et aPbOOTh . 


5.3 Non-radial heat transport 

As we mentioned in Sect. 15.11 in the case where B is nearly parallel to the surface, 
the ID approximation fails, because the heat is transported along the field lines from 
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Fig. 11 Local effective surface temperature Ts as function of the temperature Tj, at the bottom of a non- 
accreted heat blanketing envelope with Pb = 10*® g cm“^ for a neutron star with mass M = 1 AM(j, radius 
R = 12.6 km, and the dipole magnetic field with polar strength Bp = 10*® G (left panel) and 10*^ G (right 
panel). Solid lines - ID calculation with allowance for neutrino emission from the crust, dotted lines - 
neutrino emission is neglected. The upper dotted or solid curve shows at the magnetic pole, and the 
lower curve curve shows Tg at the equator. The dot-dashed curve shows the result of a full 2D calculation 
for Tg at the magnetic equator. 


the hotter surface regions outside the considered patch of the surface. Therefore the 
ID approximation overestimates the heat blanketing effect in regions with nearly tan¬ 
gential magnetic fields. For a dipole field geom etry it is an equatoria l region , whose 
width can be estimated as ^ 10% of the radius ( Potekhin et aklEoP? 1 . Since these re¬ 
gions are also the coldest ones, their contribution to the total flux is negligible. Then 
the ID approximation well reproduces the integrated observed flux. However, it is not 
the case for magnetars, which may have a complex field geometry. A 2D treatment 
shows that the ID approach is rel iable in the regions where m agnetic field lines make 
a substantial angle to the surface ( Kaminker et al Jl?0 12 . 2014 ), but it predicts too low 
surface temperat ures when the tangential magnetic field dominates (see Section 3 in 


Pons et al.ll2009l) . Therefore for magnetars one must go beyond the ID approxima¬ 


tion. Complex field configurations which lack cylindrical symmetry may require the 
full 3D treatment, which has not been done yet. 


In Fig.[TT]we show results of 2D calculations in the dipole field geometry, com¬ 
pared with the ID results. In this case, we see a substantial increase of Tg at the 
magnetic equator. This effect is especially pronounced for the superstrong field on 
the right panel. In AppendixlBTlwe give an analytical approximation to the Tj, - re¬ 
lation in the case of a strong magnetic field, including the effects of neutrino emission 
from the crust. 
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6 Thermal evolution of neutron stars 


6.1 Cooling scenarios 


Several tens of seconds after birth, the protoneutron-neutron star has lost its excess 
lepton content, i t has finished its residual con t raction and becomes transparent to neu- 


lepton content, i t has tinished its residual con t raction and becomes transparent to neu¬ 
trino emission ( Burrows and Lattimer 1986t Pons et al. 19991 Roberts! 2012 ). Soon 


after that, the temperature distribution in the highly conductive stellar core reaches 
equilibrium, which is preserved thereafter throughout the star lifetime (except dur¬ 
ing short periods after catastrophic phase transitions in the core postulated by certain 
hypothetical models). 

In the initial cooling stages, the stellar crust is hotter than the core, which is 
rapidly cooled down by the copious neutrino emission. The cooling wave reaches the 
surface within 10-100 years; thereafter, the star cools down in the quasistationary 
regime. Since all currently observed neutron stars are at least several centuries old, 
they should be in the state of quasistationary cooling, except during transient events 
with signihcant energy release in the crust or the ocean discussed below. 

Cooling in the quasistationary regime goes through two major stages. The first, 
neutrino cooling stage lasts ^ 10^ years. During this period, the core cools mostly 
via neutrino emission. The second, photon cooling stage begins when the low tem¬ 
perature of the core makes the neutrino energy loss es smaller than the losses du e to 
electromagnetic radiation from the surface (see, e.g., Yakovlev and Pethickl2004 , and 
references therein). This occurs at the age of r; 10^ years, depending on the particular 
stellar model and local conditions. 

A theoretical cooling curve of an isolated neutron star, which shows the photon 
luminosity of the star Lph or its temperature as a function of age f, depends on the 
stellar mass M, on the model of superdense matter in the core, which in particular, 
determines the intensity of neutrino emission and the EoS (and hence the stellar radius 
R), and on the properties of the envelopes. The latter include the thermal conductivity, 
which determines Lph at a given internal stellar temperature, the neutrino luminosity 
Qv in the stellar crust, and the intensity of heating sources H. For highly magnetized 
neutron stars, the cooling curve also depends on the magnetic field B (on both its 
strength and configuration), since it affects the microphysics (conductivities, EoS, 
specihc heat, etc.). Therefore, in general, the thermal evolution equations (|4|i should 
be supplemented by the equations that describe evolution of the magnetic field and 
electric currents i n the star, which leads to the thermomagnetic evolution scenarios 
(see the review by Mereghetti et al.ll2015 and references therein). 

By comparing theoretical cooling curves with the observed Lph and t of isolated 
neutron stars, one can eventually place bounds on the theoretical models of super- 
dense matter. At contrast, most neutron stars in binary systems have an additional 
source of energy (accretion) and an additional source of X-ray radiation (accretion 
disk), often much more powerful than the surface thermal emission. For this reason, 
they cannot be used to test cooling models. 

The theoretical cooling scenarios are currently divided into two main classes: 
“minimal cooling” and “enhanced cooling”. The enhanced cooling implies fast neu¬ 
trino emission processes, such as Durca reactions (Sect. 12.6.21 i. whereas the min- 
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imal cooling does not include such processes, but may be enhanced at the epoch 
of the onset of the baryon superfluidity via the PBF neutrino emission mechanism 
(Sect. I3.31 l. which helps t o explain the variety of the observed s urface temperatures 
of cooling neutron stars ( Gusakov et al.l 2004 : Page et al. 20041) . A spectacular ex¬ 
ample is the neutron star CXO J232327.9H-584842 in the Cassiopeia A supernova 
remnant, dubbed Cas A NS, w hich shows an unexpectedly appreciable tem pera¬ 
tur e decline during sev eral years ( Heinke and Holl^lOl: Elshamoutv et al.l[2013 ) (but 


see 


Posselt et al.l 1201 3l for tentative alternative interpretations of t he observations) . 


This decline can be comfort ably explained b y the PBF emission ( Page et ^ 2011 


Shternin et al. 2011 ; see also Ho et al. 20151 for a recent analysis including modern 


observational data). 


6.2 Heating mechanisms 


In the course of their evolution, some regions of neutron stars may not only cool 
but also heat up by different mechanisms. For instance, the polar cap of a pulsar can 
be heated by a stream of electrons or positrons moving along open field lines from 
the magnetosphere. The temperature of “hot spots” produced by this additional heat 
deposited onto t he stellar surface may be mu ch higher than the average temperature 
of the star (e.g., Greenstein and Hartkdll983|). Non-uniform heating processes occur 
also during accretion episodes (e.g.. ilnogamov and Sunvaevll20ld) . The hot polar 
caps emit much more intense X-rays than the remaining surface; as a result, such 
neutron stars become X-ray pulsars. Pulsed X-ray radiation is also observed from 
thermonuclear expl osions of accreted matter at the surface of a rotating neutron star 
(see, e.g., review by Strohmaver and Bildstenl[2()()^ . 

On the other hand, a neutron star may also b e heated from inside, for example due 
to dissipation of a strong magnetic field (e.g ., Mrralles. Urpin and Konenkov! 1998t 
Urpin and Konenkov 2008 : Pons et aP 2009 ). It has been suggested that dissipation 
of superstron g magnetic fields m ay be responsible for the h igh effective temperatures 
of magnetars ( ThompsonlbOOI ; see Mereghetti et al.ll2015l for a revie w). Noticeable 


liberation of energy in the crust may also oc cur during starquakes (iHaensel et al 


1990l: Franco et al. 2000; Perna and PonsIfioTll) . Two heating processes related to the 
secular spin-down of the star have also been proposed; vorte x creep, the dissip ative 
motion of superfluid vortices through the neutron star crust (lAlpar et alJI 19841) . and 
rotochemical heating, the energy released by non-equilibrium beta deca ys due to the 
slow c ontraction of the neutron star as i ts centrifugal force decreases (iReisenegger 
1995 ). Gonzalez and Reiseneggei ( 2010l) performed a comparative study of several 
heating mechanisms and found that the rotochemical heating and vortex creep can be 
most important for classical and millisecond pulsars. Both processes, albeit model- 
dependent, can keep millisecond pulsars at a surface temperature Teff ^ 10^ K. 

Another class of neutron stars undergoing heating episodes are quasipermanent 
transients, i.e., those soft X-ray transients (SXTs) whose active and quiescent periods 
last a few years or longer. During high-state accretion episodes, compression of the 
crust under the weight of newly accrete d matter results in deep crustal he ating, driven 
by exothermic nuclear transformations ( Haensel and Zdunikll990[ 2008 ). For a given 
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neutron star model, one can calculate the heating curve, that is the dependence of the 
equilibrium accretion-free on the accretion rate averaged over a large period of 
time. There is a close correspondence between the theory of therma l states of tran 


sientl y accreting neutron stars and the theory of neutron star cooling (I Yakovlev et al 


20031). Comparing the heating curves with a measured equilibri um Tpff value, one can 
constrain parameters connec t ed to properties of dense matter (lYakovlev et al.ll2004 


constrain parameters connec t ed to properties ot dense matter YaKovlev et al.N/UU4 ; 
Levenfish and Haensel 2007t Ho 120111) . Wiinands et alJ ( 2013h discuss prospects of 


application of such analysis to various classes of X-ray transients. The SXTs that 
have recently turned into quiescence a llow one to probe th e state of the neutron-star 
crust with the observed decline of 7i.ff. lBrown et al.l (119981) suggested that during this 
decline the radiation is fed by the heat that was deposited in the crust in the preced¬ 
ing active period. Such cooling is independent of the details of the star structure and 
composition and therefore its analysis directly yields information on the physics of 
the crust. Observations of several sources can be interpreted in terms of this hypoth¬ 
esis and yield cons traints to the heat conductivity in the neutron-star crust, as, e.g., 
for KS 1731-260 dShternin et alJl2007t iBrown and Cummingll2009l), XTE J1701- 


462 (Fridriksson et alJl201 1 : Page and Reddvll20 13 ). EXO 0748-676 ( Turlione et al 

201-51) . The conductivity proves to be rather high, which means that the crust has a 
regular crystalline (not amorphous) structure. On the other hand, there are similar 
objects which display variations of thermal flux that do not conform to the thermal- 
relaxation scenari o, which may be caused by a residual slow accretion on the neutro n 
star in quiescence ( Rutledge et alJE002 : Coti Zelati et al.lEoi4 : Turlione et alJl2015 ). 


6.3 Thermal luminosities of isolated neutron stars 


The inferred effective temperature depends on the choice of the emission model 
(blackbody vs. atmosphere models, composition, condensed surface, etc. - see Pa¬ 
per I), which typically results in variation of Tj-ff by a factor si 2 — 3, and it has even 
larger theoretical uncertainties in the case of strong magnetic fields. In addition, pho¬ 
toelectric absorption in the interstellar medium further constitutes a source of etTor 
in temperature measurements, since the value of the hydrogen column density Nn is 
correlated to the temperature value obtained in spectral fits. Different choices for the 
absorption model and the metal abundances can also yield different results for the 
temperature. Last, in the case of data with few photons and/or strong absorption fea¬ 
tures, the temperature is poorly constrained by the fit, adding a large statistical error 
to the systematic one. 

Because of these uncertainties, the luminosity may often be a better choice to 
compare data and theoretical cooling models. Since it is an integrated quantity, it av¬ 
erages effects of anisotropy and the choice of spectral model. The main uncertainty 
on the luminosity is often due to the poorly known distance to the source. In many 
cases, the distance is known within an error of a few, resulting in up to one order 
of magnitude of uncertainty in the luminosity. In addition, the interstellar absorption 
acts predominantly in the energy band in which most of the middle age neutron stars 
emit (£ < 1 keV). Clearly, hottes t (magnetars) or closest (XINSs) sources are easier 
to detect (see Vigano et al. 20131 and Paper I). Similarly to the case of the tempera- 























































46 


A.Y. Potekhin, J.A. Pons, D. Page 


Table 3 Cooling neutron stars, is the characteristic age, is the kinematic age, and /x is the unabsorbed 
flux in the 1-10 keV band. The range of luminosities L includes both statistical and distance errors; for 
strongly absorbed sources (i.e., most magnetars) a minimum arbitrary factor of 50% uncertainty is assumed 
to account for systematical model-deoendent uncertainties. Data have been taken fromlVigano et alJi2013h 
(see references therein and the online catalog in http: //www.neutronstarcooling, inf o.) 

Source 

logioClc) 

[yr] 

log 10 (Ik) 

[yr] 

logio(/x) 

[erg cm“^ s“^] 

d 

[kpc] 

log 10 (i-) 

[erg/s] 

CXOU J185238.6+004020 

8.3 

3.7-3.9 

-12.3 

7.1 

33.5-33.7 

IE 1207.4^5209 

8.5 

3.4-4.3 

-11.8 

0.8 

33.0-34.0 

RX 10822^300 

8.3 

3.5-3.6 

-11.3 

2.2±0.3 

33.5-33.7 

CXO J232327.9+584842 

- 

2.5 

-11.8 

3 4+^-3 

0.1 

33.4-33.6 

PSR J0538+2817 

5.8 

»4.6 

-12.1 

1.3±0.2 

32.7-32.9 

PSRB1055-52 

5.7 

- 

-13.4 

0.73±0.15 

32.2-32.6 

PSR J0633+I746 

5.5 

- 

-12.5 

^•^'^-0.08 

31.6-32.5 

PSR 61706^4 

4.2 

- 

-12.1 

0.6 

31.7-32.1 

PSR 60833^5 

4.1 

3.7-4.2 

-10.5 

0.28±0.02 

32.1-32.3 

PSR B0656+14 

5.0 

«4.9 

-12.6 

0.28±0.03 

32.7-32.8 

PSR B2334+61 

4.6 

«4.0 

-14.0 

J.i 24 

30.7-32.1 

PSRJ1740+1000 

3.1 

- 

-13.8 

1.4 

32.1-32.2 

PSRJ1741-2054 

5.6 

- 

-12.5 

0.8 

30.4-31.4 

PSR J0726-2612 

5.3 

- 

-14.0 

1.0 

32.1-32.5 

PSR J1119-6127 

3.2 

3.6-3.9 

-13.0 

8.4±0.4 

33.1-33.4 

PSRJ1819-1458 

5.1 

- 

-12.6 

3.6 

33.6-33.9 

PSRJ1718-3718 

4.5 

- 

-13.2 

0.0 

32.8-33.5 

RX J0420.0-5022 

6.3 

- 

-17.8 

0.34 

30.9-31.0 

RXJ1856.5-3754 

6.6 

5.5-5.7 

-14.4 

0.12±0.01 

31.5-31.7 

RXJ2143.0+0654 

6.6 

- 

-13.1 

0.43 

31.8-31.9 

RXJ0720.4-3125 

6.3 

5.8-6.0 

-13.3 

0.02 

32.2-32.4 

RXJ0806.4-4123 

6.5 

- 

-13.4 

0.25 

31.2-31.4 

RXJ1308.6+2127 

6.2 

5.9-6.1 

-12.1 

0.50 

32.1-32.2 

RXJ1605.3+3249 

- 

5.7-6.7 

-13.0 

0.35±0.05 

30.9-31.0 

IE 2259+586 

5.4 

4.0-4.3 

-10.3 

3.2±0.2 

35.0-35.4 

4U 0142+614 

4.8 

- 

-9.8 

3.6±0.5 

35.4-35.8 

CXO 1164710.2^55216 

5.2 

- 

-12.2 

4 

^•^-1.0 

33.1-33.6 

XTE 11810-197 

4.1 

- 

-11.7 

3.6±0.5 

34.0-34.4 

IE 1547.0-5408 

2.8 

- 

-11.5 

4.5±0.5 

34.3-34.7 

IE 1048.1-5937 

3.7 

- 

-10.8 

2.7±1.0 

33.8-34.5 

CXOU 1010043.1-721 

3.8 

- 

-12.5 

60.6±3.8 

35.2-35.5 

IRXS 1170849.0^00910 

4.0 

- 

-10.4 

3.8±0.5 

34.8-35.1 

CXOU 1171405.7-381031 

3.0 

«3.7 

-11.4 

13.2±0.2 

34.9-35.2 

IE 1841-045 

3.7 

2.7-3.0 

-10.4 

9.6+«-6 

35.2-35.5 

SGR 0501+4516 

4.2 

w4 

-11.3 

0.5 

33.2-34.0 

SGR 1627^1 

3.3 

«3.7 

-11.6 

11.0±0.2 

34.4-34.8 

SGR 0526-66 

3.5 

«3.7 

-12.0 

49.7±1.5 

35.4-35.8 

SGR 1900+14 

3.0 

3.6-3.9 

-11.1 

12.5±1.7 

35.0-35.4 

SGR 1806-20 

2.6 

2.8-3.0 

-10.6 

13.0^t!! 

35.1-35.5 

SGR 0418+5729* 

7.6 

- 

-14.0 

2.0 

30.7-31.1 

Swift 11822.3-1606* 

6.2 

- 

-11.5 

1.6±0.3 

32.9-33.2 


Notes. *The source has been recently discovered in outburst and it could have not yet reached the quies¬ 
cence level. 
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ture, the choice of different models of absorption and chemical abundances can yield 
additional systematic errors on the luminosity. However, for the worst cases, the rela¬ 
tive error is about 30%, making it usually a secondary source of error compared with 
the distance. 

In Table [3 we summarize the properties of cooling neutron stars@ and in 
Fig. \T2\ we c ompa re the observational data to theoretical cooling curves, from 


Vigano et al.l (12013 1. Here, the theoretical results are computed by a finite differ- 


ence method for 2D (axisym metric) stellar configurations, usin g the SLy EoS mo del 
( Douchin and Haensell 200 ih at p > Pd^p and the BPS EoS ( Bavm et al. 1971 ) at 
p < Pdrip. The high Durca threshold of the SLy EoS has been artihcially lowered for 
illustrative purpose to p = 10*^ g cm ~^, corresponding to the central density of a star 
with M =l A2 Mq (see Viganol2013 for d etails). Eor superfluid gap energies, the phe¬ 
nomenological model of lHo etal.N2012h has been adopted. The other microphysics 
input is the same as in Sects. EPS] 

In the upper panel of Eig.[T2]we show cooling curves for non-magnetized neutron 
stars with masses ranging between 1.10 and 1.76 Mq (lines from top to bottom). 
After Ri 100 yr, low mass stars (M < IAMq) are brighter than high mass stars. Eor 
the high-mass family, M >l AMq, the Durca processes in the central part of the star 
result in fast cooling before one hundred years. Within the low-mass family, cooling 
curves are similar at early ages (< 100 yr). The differences at f ~ 10^ — 10^ yr are due 
to the delayed transition of neutrons in the core to a superfluid state, which activates 
the PBE neutrino emission. After the effect of the transition to a superfluid core is 
finished, at f > 10^ yr, cooling curves for low-mass neutron stars tend to converge 
again, following the same curve independently of the mass. 

We see that luminosities of some objects in the upper panel of Pig.[T2]are system¬ 
atically above the theoretical cooling curves. Eor the CCOs this discrepancy can be 
eliminated by considering accreted (more heat-transparent) blanketing envelopes, as 
the lowest dashed line in the lower panel of Eig.[T2]demonstrates. However, the high- 
B objects still remain systematically hotter than what the theory can explain at Z? = 0. 
This provides strong evidence in favor of the scenario in which magnetic field de¬ 
cay powers their larger luminosity. In the lower panel we compare the observational 
data to theoretical cooling curves for different values of the initial magnetic field up 
to 3 X 10*^ G. The most relevant effect of the inclusion of the magnetic held is that 
it allows to explain objects with high luminosities. Magnetic helds B > lO*"*^ G are 
strong enough to noticeably heat up the crust and power the observed X-ray radiation. 
Another important difference is that the cooling timescale for strongly magnetized 
objects is several times larger than for the weakly magnetized neutron stars. 


7 Conclusions 

We have considered the basic physical ingredients needed for theoretical modeling 
neutron-star thermal evolution and briehy reviewed some recent results on cooling 
of magnetized neutron stars. The physics behind such thermal evolution is extremely 

^ A regularly updated online catalog can be found at http: //www.neutronstar cooling, inf o, with 
abundant links to references for each source. 
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Fig. 12 Comparison between observational data and theoretical cooling curves (from lVigano et alJ2013l) . 
The obs ervational estimates of (errorbars) or constraints on (arrows) the age and thermal luminosity corre¬ 
spond to lVigano et alJ j2Q13h and Table[3] The abbreviations in the legend mark different classes of neutron 
stars with measured thermal radiation (MAG - magnetai* candidates, XINS - X-ray isolated neutron stars, 
HB - high-5 radio pulsars, RPP - rotation powered pulsars, CCO - central compact objects; see Paper I). 
Upper panel: non-magnetic neutron stars with iron envelopes, withM = (1.10, 1.25, 1.32, 1.40, 1.48, 1.60, 
1.70, 1.76) Mq (lines from top to bottom). Lower panel: a neutron star with M — \AMq and R— 11.6 
km, and three different cases with initial mag netic field at the pole B - 0, 3 X 10^4 G, and 3 x 10^^ G. The 
magnetic field topology is that of Model A in IVigano et alj j2013h (crustal confined). We show results for 
iron envelopes (solid lines) and hydrogen envelopes (dashed lines). 
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rich. Clearly, we could not consider it in depth in a single review paper. However, 
the information that we have given, together with the references to the formulae and 
online resources elsewhere, should be sufficient to build a neutron-star cooling model 
involving only the simplest assumptions. We considered the basic equations that gov¬ 
ern the mechanical and thermal structure of a neutron star and its thermal evolution, 
the main contributions to the physical quantities that enter these equations - namely, 
EoS and heat capacity, thermal conductivity, neutrino emissivity, the effects of baryon 
superfluidity and proton superconductivity and of strong magnetic fields. In addition, 
we present a novel fit to the relation between the internal and external temperatures 
and heat fluxes in the blanketing envelope, which includes the effects of neutrino 
emission from the crust and the effects of non-radial heat transport. 

In this paper we have restricted ourselves by the npejJ. matter, without either 
hyperons or “exotic” models that involve hyperon condensates, quark phases, mixed 
phases, or phase transitions. We hope that an interested reader should be able to study 
these issues in depth following the literature references that we have provided. We 
have not considered also the equations of magnetic-field evolution, coupled to the 
thermal evolution, which is espec ially important in magn etars. These equations are 
given, for instance, in the paper by Mereghetti et alJ ( 2015h in this volume, where ori¬ 
gin, evolution, and observational manifestations of magnetars are reviewed in depth. 
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Appendices 

A. Electron thermal conductivities 

In this Appendix, we briefly overview the physics of electron heat conduction in the neutron-star envelopes, 
which is the most important heat conduction mechanism as regards the neutron-star thermal evolution, in 
the case of B = 0. The magnetic field effects on the heat conduction are considered in Sect. EM 


A. 1 Weakly degenerate electron gas 

In the case o f non-degenerate and non-relativistic electrons jSpitzer and Harmlll953l: [Braginskij|l958t 
ISpitzeilll962h . the effective energy-averaged electron-ion collision frequency is 

(A.l) 

where A^i is the Coulomb logarithm. In the considered case A^i is a slowly varying function of density 
and temperature. Its precise value depends on the approximations used to solve the Boltzmann equation, 
but its order of magnitude is given by the elementary theory, where the Coulomb collision integral is 
truncated at small and large impact parameters of the electrons. Then A^i ~ ln(rniax/^min)i where r^ax and 
rmin are the maximum and minimum electron impact parameters. The parameter rmax can be set equal 
to the Debye screening length, r~l^ — 47i{ne -\~Z^n[)e^/T. The second parameter can be estimated as 


4 [271 ^ 
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= max(Ae, rd), where (defined in Sect. 12.3t limits rmm in the high-temperature regime (where 
the Bom approximation holds), and — Ze^/T is the classical closest-approach distance of a thermal 
electron, which limits rmin in the low-temperature, quasiclassical regime. 

A similar effective frequency 




8 fir e" 
3V^ 7^”' ' 


(A.2) 


characterizes the efficiency of the ee collisions. If Aee ~ A^i, then v^JVee ~ Z, therefore for large Z the ei 
collisions are much more efficient than the ee collisions. 


A.2 Strongly degenerate electron gas 
A.2.7 Electron-ion scattering 

The thermal conductivity of strongly degenerate electrons in a fully ionized plasma is given by Eq. i\2\ 
with a — /I. In order to determine the effective collision frequency that enters this equation, we use the 

Matthiessen rule v = v^i + Vee- 

_ The effective electron-ion collision frequency can be written in the form jLedl95(]l : lYakovlev and UrpinI 

119801) 

^Zm*e^Ae\ ZAp\ \/\-\-x} 

Vei =-. (A.3) 

5.7xl0-^^s 

iLed 1195Q|) gave an estimate of the Coulomb logarithm = InCrma x/^mm), with the minimum impact 
parameter rmin — fi/lp-p and the maximum impact parameter r^ax = ^;. lYakovlev and UrpinI<198(1) calcu¬ 
lated the conductivities for relativistic degenerate electrons, neglecting electron screening, and obtained a 
more accurate estimate rmax ~ OAai in the liquid regime. In the so lid regime, where the electrons scatter 
on phonons (collective ion excitations). lYakovlev and UrpinI jl98Ql) obtained different approximations for 
the two distinct cases, 6 )d <T <T^ and T < 6)d- 

IPotekhin et al.l j 19991) derived a unified treatment of the electron conductivities in the Coulomb liquid 
and solid and described both regimes by Eq. i|a 3). Then qualitatively, by order of magnitude, A^j ~ 1 
in the ion liquid, and Ae\ ~ T/Tm in the Coulomb solid with a melting temperature T^. The effects of 
multiphonon scattering, electron screening, and non-Bom coiTections, have been taken into account, and 
the Coulomb logarithms in both liquid and solid phases have been fitted by a single analytical formula. A 
Fortran code and a table of thermal conductivities, based on this formalism, are available online[3 

At the conditions typical for the envelopes of neutron stars, the electron-phonon scattering pro¬ 
ceeds mainly via the Umklapp processes, where t he wave vector con'espondin g to the change of elec¬ 
tron momentum lies outside the first Brillouin zone. lRaikh and YakovlevI jl982l) noticed that if T < Ty = 
TpZ^/^af x/1 +jr^/3jCr, then the Umkla pp processes occur less often (“freeze out”). Then the scattering rate 
decreases. iRaikh and YakovlevI jl982h assumed an extremely strong (exponential) decrease. This implied 
that at T < 7u the conductivity would be in practice d etermined by impur ities and s tructure defects of 
the lattice, rather than by the electron-phonon scattering jGnedin et alj200lh . However. IChugun^ j2Q12h 
showed that distortion of electron wave functions due to interaction with the Coulomb lattice destroys this 
picture and strongly slows down the increase of the conductivity. As a result, the conductivities in neutron 
star envelopes can be treated neglecting the “freezing-out” of the Umklapp processes. 

A.2.2 Electron-electron scattering 

Although the electron-ion scattering is usually most important for degenerate pla smas, the ele ctron- 
electron scattering still can be non-negligible for relatively light elements (Z < 10) jLampdll968h. The 
expres s ion of Vg^ for the relativist ic degenerate electrons at T -C Tp was obtained bv iFlowers and It5] 
jl976l).IU rpin and Yakovle^ j 19801) extended it to higher temperatures, where TJ, < 7 fip. 

IShtemin and YakovlevI j2QQ6h reconsidered the problem including the Landau damping of transverse 
plasmons, neglected by the previous authors. This effect is due to the difference of the components of 
the polarizability tensor, responsible for screening the charge-charge and current-current interactions: 
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http://www.ioffe.ru/astro/conduct/ 
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the transverse current-current interactions undergo “dynamical screening.” IShtemin and Yakovlev! j2Q06h 
showed that the Landau damping of transverse plasmons strongly increases v^e in the domain of Xp > 1 
and r Tp and presented a new fit to (also implemented in the code referenced in footnote[4). 


A.3 The case of intermediate degeneracy 


In the case where the electron gas is partially degenerate, that is T ~ the thermal and electrical con¬ 
ductivities determined by the electron-ion scattering are satisfactorily evaluated by the thermal averaging 
procedure [Eq. {33) in Sect. 1332]. For conductivities determined by the elecb'on-electron collisions, there 
is no such averaging procedure, but we can use an interpolation between the two limiting cases, 


V -yde. l+625(r/eF)^ 

" ‘‘ l+25r/eF + 271(r/eF)5/2’ 

A satisfactory accuracy of this interpolation has been verified bv ICassisi et alJi2007ll . 


(A.4) 


A.4 Impurities and mixtures 


If the plasma in an envelope is not a pure substance of a single chemical element, then the effective collision 
frequency Vgj should be modified. The required modificati on can be different, dep e nding on the state of 
the pla sma and on the amount of imp urities. For example. iFlowers and Itohl jl976l) . [Yakovlev and Umii] 
jl98Ql) . and lltoh and Kohvamal jl993h considered electron scattering by charged impurities in a Coulomb 
crystal. If the fraction of impurities is small and they are randomly distributed, then electron-impurity 
scattering can be treated as scattering by charge fluctuations, controlled by the impurity parameter Q — 
((Z— (Z))^), where (Z) = Yj — nj/Y^jiij is the number fraction of ions of the yth kind, and 

Zj is their charge number. Then, using the Matthiessen rule, one can obtain Vgi as a sum of the terms 
corresponding to the electron-phonon scattering in a homogeneous lattice and to the electron scattering 
by charge fluctuations. The effective relaxation time for the latter term is given by Eq. <A.3t with ZAgi 
replaced by where the Coulomb logarithm Aj depends generally on j. Neglecting 

the differences between the Coulomb logarithms, one can thus simply replace Z by Ql{Z) in Eq. iA3) to 
estimate the conductivity due to electron scattering by chai'ged impurities. 

An alternative approach is relevant when there is no dominant ion species which forms a crystal (e.g., 
in a liquid, a gas, or a glassy alloy). In this case, one can use Eq. {a 2) with Z'^jAgi replaced by Y. i Z^n jA . 
An approximation to Aj based on the plasma “additivity rule” has been suggested bv IPotel^in et alJfl999h . 
Neglecting the differences between the Coulomb logaiithms, one anives at Eq. {a 3) with Z replaced by 
yJ{Z-). If tabulated conductivities Kj for pure substances are used, then the best agreement with calcula¬ 
tions based on the “additivity rule” is usually given by the estimate 


_ {kz) 
LjYjZj - (Z) • 


(A.5) 


B. Temperature relations for envelopes of neutron stars with magnetic fields 

Here we present an analytical fit to the temperature distribution over a surface of a neutron star with 
a non-accret ed envelope and a dipole magnetic field. We have chosen Ph — 1 0*^ g cm“^ and used the 
BSk21 EoS JPearson et alJ201^ in th e parametrize d form ^Potekhin et alJ20T3) . The numerical data have 
been produced with the 2D code of IVigano et alj (203) for 5 values of internal temperature T\j from 
10^ K to 10^ K, 5 values of the magnetic field at the pole Bp from 10‘^ G to 10‘^ G, and 20 values of 
magnetic colatitude 9 at the surface of the neutron star from 0 to ;r/2. The use of the 2D code corrects 
the temperature distribution near the magnetic equator, because the non-radial heat flow increases the 
equatorial as compared to the ID model that was employed previously (see Fig. llll in Sect. 15.2t . These 
data ha ve been supplemented with more detailed calculations at the magnetic pole (6 = 0) using the ID 
code of IPotekhin et all l2007^ for 36 values of T), from 10®'^ K to 10*° K and 9 values of Bp from 10* * G 


























52 


A.Y. Potekhin, J.A. Pons, D. Page 


to 10^^ G. An important difference from the old results is the inclusion of the neutrino emission from the 
crust, which is especially important for the magnetars (see Sect. [5). Bec ause of the 2D treatm ent and the 
allowance for neutrino emission, the new fit supersedes the previous one jPotekhin et alJ2QQ5) . whenever 
B > 10^^ G or 7b > 10^ K. We stress that its use is restricted by non-accreted (i.e., composed of heavy 
chemical elements) envelopes in the range of 10^'^ K < 7b < 10*^ K and < 10*^ G, which is covered 
by the underlying numerical data. For envelopes with B < 10^^ G (either non-accreted or accreted), the 
previous fit can be used, however the surface temperature (but not the flux at the inner boundary, Fb — 
see item 4 below) should be limited for hot stars according to Eq. El below. 

The fit consists of 3 stages: (1) an expression for the surface temperature at the magnetic pole, Tp, as 
function of 7b, g, and B^\ (2) an expression for the ratio of the polar to the equatorial surface temperatures, 
T’p/Teq; (3) an expression for the dependence of 7’s on the magnetic colatitude 9. Since the thermal conduc¬ 
tivities for quantizing magnetic fields (Sect. |4A2t are known for the electron-ion but not electron-electron 
collision mechanism, we multiplied by a correction factor, obtained numerically from a comparison of 
the results of thermal-structure calculations with and without the ee collisions at F = 0. At the end of this 
Appendix we suggest a recipe for relating the flux Fb at the bottom of the heat-blanketing envelope to 
temperature Fg and thereby to Fb. 

1 . At the magnetic pole, the effective surface temperature, neglecting neutrino emission from the crust, 
is approximately given by the expression 

4 °^ = [gl4(7’lV(l+0.15^/^)^o^)]‘^'‘ X 10« K, (B.I) 

where 

^0 = (15.77’5^''^ + 1.36^9)°■”'’^ = a = 0.337/(1+0.02^^), (B.2) 

Fg — Fb/IO^ K, and B \2 — Sp/10*^ G. The limiting temperature, at which Fp(Fb) levels off due to the 
neutrino emission from the cnast is approximately given by 

= (5.2gOf +0.093X 10« K. (B.3) 

The corrected surface temperature at the pole, which takes this limit into account, is reproduced by the 
expression 

[l + (B.4) 

2 . The ratio of the polar to equatorial surface temperatures can be roughly evaluated as 


(1230r9p^fii2^1+2flf2 0.0066Bff 

Tcq ^ (Bi2 + 450r9 + 119Bi2 79)'^ ^ r,'^^ + 0.00258Bf^^ 


The numerically calculated Fp/Feq ratio has a complex dependence on Fb and B at B > 10^^ G. In order 
to keep our fitting formulae relatively simple, we do not reproduce these oscillations, but instead force the 
ratio ED to converge to some average value atB^ 10^^ G. The numerical data oscillate in a complicated 
manner around this average, with deviations reaching up to 35%. For smaller fields, F < 3 x 10^^ G, 
Eq. <B3) reproduces the numerical data with typical en'ors of several percent (up to 10%). Note that 
these significant deviations affect only nearly tangential field case, viz. the equatorial region, which is 
substantially colder than the rest of the surface. Therefore its contribution to the observed flux is usually 
not very important. 

3 . Finally, the dependence of the surface temperature on the magnetic colatitude 9 is approximately 
described by the expression 


Ue)-T,^ 

Tp ~ Teq 


(1 + +fl2)cOS^0 
1+^1 COS 0+^2 COS^ 9 ' 


where ai — 



02 


10Fi2 

+0.1 


(B.6) 


4 . Note that the outer boundary condition to the thermal evolution equations 0 involves the relation 
between the heat flux density Fb through the boundary at p = Pb and the temperature Fb at this boundary. 
In the absence of the neutrino emission from the crust, this boundary condition is directly provided by 
the Fb-Fs relation, because in this case (in the plane-parallel approximation) Fb = CJshT^- It is not so if a 
significant part of the energy is carried from the outer crust by neutrinos. In this case we suggest to evaluate 
the flux through the boundary by the relation Fb = <7 sbF^^, where F* is given by the above approximations 
for Fg, but without the correction ED- 
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